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ABSTRACT 

/ 

This  study  attempts  to  evaluate  the  influence  of 
soil  conservation  practices  on  the  surface  hydrology  in 
the  Galena  watershed,  from  1940  through  1987.  Poor  soil 
conservation  practices  within  the  Galena  watershed  from 
the  mid-1800s  until  the  mid-1900s  increased  surface 
runoff  while  decreasing  baseflow.  The  poor  soil 
conservation  practices /inf luenced  soil  erosion,  flood 
magnitudes,  channel  morphology,  and  sedimentation 
processes. 

The  influence  of  poor  soil  conservation  practices 
within  the  Galena  watershed  reached  a  maximum  during  the 
early  1900s.  Efforts  begun  by  the  Soil  Erosion  Service 
in  the  mid-1930s  began  to  reduce  surface  runoff  by 
educating  farmers  in  proper  soil  conservation  methods. 
Aerial  photography  analysis  indicates  a  steady  increase 
in  strip  cropping  and  contour  farming  since  1937.  The 
increased  use  of  soil  conservation  practices  within  the 
Galena  watershed  have  decreased  surface  runoff  more  than 
20%  and  increased  the  baseflow. 

The  influence  of  soil  conservation  was  evaluated  by 
controlling  climatic  variables  through  an  analysis  of 
similar  storm  events.  Statistical  analyses  indicates' 
similarity  among  the  storm  events,  reducing  the  climatic 
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influence.  Hydrographs  corresponding  to  the  selected 
storms  were  analyzed  and  an  estimation  of  surface  runoff 
and  baseflow  calculated  for  each  storms  throughout  the 
study  period.  Overall  results  indicate  a  general 
decreasing  trend  in  surface  runoff  with  a  corresponding 
increasing  trend  in  baseflow  during  the  study  period. 


The  primary  objective  of  thie  study  is  to  compare 
storm  runoff  hydrology  as  it  changed  f roe  the  early 
1940s  through  the  late  1980s.  A  specific  attempt  is 
made  to  evaluate  the  influence  that  soil  conservation 
practices  have  had  on  the  surface  runoff  hydrology  of 
the  intensively  cultivated,  topographically  steep  Galena 
watershed  upstream  of  the  U.S.  Geological  Survey 
(U.S.G.S.)  stream  gage  (105415000)  at  Buncombe, 

Wisconsin  from  1940  through  1967.  Specific  goals  of 
this  study  are  to  quantitatively  determine  the  influence 
that  cultivation  and  conservation  practices  have  had  on 
the  following: 

1.  the  surface  runoff  fraction  of  total 
runoff; 

2.  the  base  flow  fraction  of  total  runoff; 

and, 

3.  the  lag  time  between  precipitation  events  and 
peak  river  discharge. 

The  study  area  is  a  medium  sized,  intensely 
cultivated,  unglaciated  watershed  located  in 
southwestern  Wisconsin  (see  Figure  1) .  The  agricultural 
land  use  in  this  watershed  changed  significantly  from 
the  early  1940s  to  the  late  1980s  due  to  fanners 
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Figure  1.  The  study  area  is  that  portion  of  the 
watershed  upstream  of  the  Buncombe  gaging  station, 
small  portion  lies  within  Illinois. 
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adopting  soil  conservation  practices. 

Migration  into  the  Midwest  U.S.  during  the  19th  and 
early  20th  centuries  may  be  the  most  flagrant  example  of 
land  abuse  by  poor  cultivation  practices  (Butzer, 

1974) .  Reckless  cultivation  practices  emphasized 
cultivating  all  available  slopes,  plowing  perpendicular 
to  slope  contours,  and  stressing  monocultures  with 
practically  no  intercropping.  In  approximately  150 
years  the  agricultural  soil  resources  within  the  Midwest 
were  reduced  by  perhaps  50%  (Butzer,  1974)  .  During  the 
1930s  in  the  Midwest  land  abuse  began  to  decline  through 
educational  efforts  by  the  U.S.  Soil  Erosion  Service  and 
implementation  of  agricultural  conservation  practices. 

The  first  intensive  white  settlement  of  southwest¬ 
ern  Wisconsin  began  in  the  early  19th  century  with 
the  discovery  of  rich  lead  and  zinc  deposits.  However, 
within  a  decade  or  two  it  was  the  agriculture  potential 
of  the  area  that  was  responsible  for  continued 
settlement  (Trewartha,  1940;  Blanchard,  1924) .  Mining 
was  a  dominant  activity  until  the  mid-1840s  due  to  its 
profitability,  the  uncertain  agricultural  productivity 
of  the  land,  threats  of  Indian  attacks  until  the  close 
of  the  Black  Hawk  War  in  1832,  and  federal  restrictions 
preventing  ownership  and  use  of  "mineral  lands"  for 
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agricultural  purposes  (Blanchard,  1924) .  Federal 
policies  changed  in  1446  and  agricultural  activity 
intensified  until  the  early  1900s.  Since  the  late  1890s 
agricultural  activity  in  this  region  has  remained 
intense  although  the  nature  of  land  use  has  changed 
(Trimble  and  Lund,  1982;  Blanchard,  1924). 

Many  farmers  from  western  Europe  were  accustomed  to 
gentle,  persistent  rainstorms  and  were  not  prepared  for 
the  high  intensity  rainstorms  of  the  Midwest.  Unknowing 
of  the  consequences,  farmers  cultivated  sloping  hill 
tops  and  steep  valley  side  slopes.  Such  cultivation 
worked  well  in  western  Europe  but  was  not  well  suited 
for  the  intense  rainstorms  encountered  in  the  Midwest 
and  the  easily  erodible  loess-derived  soils  of  the 
region.  As  time  progressed  and  farming  practices 
remained  unchanged,  surface  runoff  increased  and  soil 
erosion  accelerated. 

The  stock  market  crash  and  the  Dust  Bowl  drought  of 
the  early  1930s  had  a  profound  influence  on  the  Midwest 
as  well  as  the  entire  nation.  These  two  events  created 
an  awareness  of  land  resource  and  management  policies. 
Some  environmentalists  recognize  the  period  from  the 
early  1930s  until  World  War  II  as  the  Golden  Era  of 
environmental  conservation  because  of  the  many  programs 
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and  policies  then  initiated  towards  improving  the 
environment  (Petulla,  1977) .  Due  to  severe  soil  erosion 
in  southwestern  Wisconsin  caused  by  farming  practices, 
the  U.S.  Soil  Erosion  Service  (later  the  U.S.  Soil 
Conservation  Service)  and  the  University  of  Wisconsin 
established  the  first  conservation  demonstration  area  in 
1933  located  in  the  Coon  Creek  basin,  approximately  125 
kilometers  northwest  of  the  Galena  River  watershed 
(Wisconsin  Blue  Book,  1989) .  Goals  of  this  conservation 
demonstration  area  were  to  educate  farmers  in  soil 
conservation  practices  and  to  establish  demonstration 
sites  in  communities  to  disseminate  these  practices 
(Johansen,  1969) . 

The  soil  conservation  practices  taught  at  Coon 
Creek  disseminated  slowly  to  the  surrounding  areas  as 
only  233  farmers  utilized  soil  conservation  practices  in 
1939  within  a  120  kilometer  radius  of  the  Coon  Creek 
project  (Johansen,  1969) .  Aerial  photographs  of  the 
Galena  watershed  in  1937  indicate  no  cultivated  land 
under  the  common  soil  conservation  practices  promoted  by 
the  Coon  Creek  demonstration  project.  The  number  of 
Wisconsin  farmers  practicing  soil  conservation  within  a 
120  km  radius  of  Coon  Creek  increased  to  6402  by  1967 
(Johansen,  1969) .  Similarly,  aerial  photographs  taken 
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in  1985  indicate  a  substantial  fraction  of  cultivated 
land  within  the  Galena  watershed  under  soil  conservation 
practices.  Increased  utilization  of  soil  conservation 
practices  have  influenced  surface  hydrology  and  appear 
to  have  significantly  reduced  surface  runoff  in  the 
Galena  watershed  during  the  study  period. 

This  study  derives  its  primary  data  from  aerial 
photographs  and  climatic  and  discharge  records.  Aerial 
photographs  were  taken  at  sporadic  intervals  by  the  U.S. 
Department  of  Agriculture  and  U.S.G.S.  from  1937  to 
1985.  Climatic  data  from  four  weather  stations 
surrounding  the  Galena  watershed  provide  the  necessary 
precipitation  data.  The  U.S.G.S.  gaging  station  located 
on  the  Galena  River  near  the  Wisconsin  -  Illinois  border 
at  Buncombe  provided  the  necessary  discharge  data.  All 
records  for  this  study  date  from  1940. 

Agricultural  land  use  changed  significantly  from 
the  1940s  through  the  1980s.  This  study  attempts  to 
compare  the  storm  runoff  from  the  1940s  to  the  1980s 
within  the  intensely  cultivated  Galena  watershed, 
focusing  on  hydrologic  responses  to  increased  levels  of 
adoption  of  soil  conservation  practices. 
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CHAPTER  2 :  BACKGROUND 


INTRODUCTION 

Changes  in  land  use  disrupt  surface  hydrology  by 
influencing  the  runoff  processes  within  a  watershed. 

Two  principle  components  of  runoff  are  baseflow  and 
surface  runoff.  Baseflow  consists  of  groundwater  runoff 
and  delayed  subsurface  runoff  (Chow,  1964) .  The  term 
'surface  runoff'  best  describes  the  combination  of 
Horton  overland  flow,  subsurface  quick  flow,  and 
saturation  overland  flow  (Mockus,  1972) .  Horton 
overland  flow  occurs  when  rainfall  intensity  exceeds 
soil  infiltration  capacity  (Dunne  and  Leopold,  1978). 
Hydrographs  illustrate  both  baseflow  and  surface  runoff 
within  a  watershed.  Most  authors  agree  that  specific 
factors  influencing  the  runoff  process  can  be  classified 
as  either  climatic  or  physiographic  influences  (Chow, 
1964;  Wisler  and  Brater,  1963;  Dunne  and  Leopold, 

1978) .  Climatic  and/or  physiographic  changes  within  a 
watershed  often  result  in  changes  in  the  baseflow  and 
surface  runoff  components  of  a  hydrograph.  This  study 
examines  storm  runoff  hydrographs  of  the  Galena  River  at 
Buncombe,  Wisconsin  (drainage  area  323  km2)  to 
evaluate  the  influence  of  changing  agricultural  land  use 
on  runoff  hydrology  since  1940. 


The  baseflow  component  of  a  hydrograph  fluctuates 
conservatively  in  response  to  precipitation  events  in 
comparison  to  more  dramatic  fluctuations  of  the  surface 
runoff  component.  Singh  (1972)  identified  three  major 
factors  that  influence  baseflow  magnitude: 

1.  hydraulic  characteristics  of  surficial 
soils ; 

2.  groundwater  aquifer  characteristics ; 

and, 

3.  evapotranspiration  demands. 

Surficial  soil  characteristics  determine  the 
infiltration  rate  and  the  influence  of  antecedent 
precipitation  which  in  turn  determine  the  storage 
capacity  and  available  water  for  aquifers.  Groundwater 
aquifer  characteristics  affect  baseflow  by  influencing 
the  rate  that  groundwater  moves  into  stream  channels. 
Permeable  aquifers  will  transport  groundwater  more 
rapidly  into  stream  channels  than  less  permeable 
aquifers.  Evapotranspiration  demands  dictate  the 
available  water  for  baseflow.  Baseflow  usually  declines 
in  summer  due  to  increased  evapotranspiration. 

Surface  runoff  is  responsive  to  many  aspects  of 
watershed  physiographic  and  climatic  factors  (Dunne, 


1983).  Rainfall  intensity,  storm  duration,  antecedent 
soil  moisture,  soil  cover,  and  infiltration  capacity 
influence  Horton  overland  flow.  Surface  runoff  occurs 
when  rainfall  exceeds  storage  and  infiltration 
capacities.  Horton  overland  flow  occurs  frequently  in 
arid  and  semi-arid  areas  where  vegetation  is  sparse  and 
rainfalls  are  commonly  intense. 

Initial  soil  moisture,  topography,  and  the  geometry 
of  the  subsurface  flow  region  influence  the  volume  and 
rate  of  subsurface  flow.  Watersheds  characterized  by 
abundant  subsurface  flow  tend  to  produce  hydrographs 
with  low  peaks  and  relatively  flat  rising  and  recession 
limbs.  Subsurface  flow  is  dominant  in  areas  with 
permeable  soils  and  dense  vegetation. 

Saturation  overland  flow  consists  of  return  flow 
and  direct  precipitation  onto  saturated  surfaces  (Dunne 
and  Leopold,  1978) .  Large  rainstorms,  wet  antecedent 
soil  conditions,  long  hillslopes  with  gentle  gradients, 
and  shallow  soils  with  an  infiltration  capacity  that 
exceeds  the  rainfall  intensity  but  is  too  low  to  convey 
all  the  rainfall  as  subsurface  flow  influence  saturation 
overland  flow.  Areas  with  watertables  that  easily  and 
rapidly  rise  to  the  ground  surface  during  rainfall 
produce  saturation  overland  flow.  Saturation  overland 
flow  dominates  watersheds  that  are  sparsely  vegetated 
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with  thin  soils,  have  high  water  tables,  and  have  long 
gentle  concave  hillsides. 

Climatic  factors  may  be  grouped  into  three  specific 
categories  including  precipitation,  interception,  and 
evapotranspiration  (Chow,  1964) .  Precipitation  factors 
include  the  form  of  precipitation,  storm  type, 
intensity,  duration,  areal  distribution,  frequency  of 
occurrence,  direction  of  storm  movement,  antecedent 
precipitation,  and  soil  moisture.  Interception  factors 
focus  on  vegetation  composition,  age,  density  of  stands, 
and  seasonal  development.  Evapotranspiration  factors 
include  temperature,  wind,  atmospheric  pressure,  the 
nature  and  shape  of  evaporative  surfaces,  solar 
radiation,  soil  moisture,  humidity,  and  types  of 
vegetation. 

Physiographic  factors  include  physical  factors  and 
land  use  (Chow,  1964) .  Physical  factors  include 
watershed  size,  shape,  slope,  orientation,  elevation, 
stream  density,  and  related  factors  emphasizing 
infiltration  processes  and  soils.  Land  use  is  a 
separate  consideration  in  this  study  due  to  its  dominant 
hydrologic  role  during  the  study  period  of  the  1940s  - 


1980s. 


13 


RUNOFF  PROCESS 

Most  climatic  and  physiographic  factors  remain 
relatively  constant  during  this  study  period  and  have  no 
direct  influence  on  the  aspects  analyzed.  Land  use 
changes  during  this  study  period  and  has  a  significant 
influence  on  the  runoff  process. 

The  runoff  process  is  continuous  over  time  and  does 
not  have  a  beginning  or  ending  point.  Hoyt  (1942) 
describes  the  runoff  process  by  delineating  five 
phases: 

1.  the  rainless  period  just  prior  to  the 
beginning  of  rainfall  and  after  an 
extended  dry  period; 

2.  the  initial  period  of  rain  which  involves: 

a.  precipitation  in  the  channel. 

b.  interception  by  vegetation. 

c.  infiltration  into  the  soil. 

d.  temporary  retention  in  small 
depressions; 

3.  the  continuation  of  rainfall  at  variable 
rates ; 

4.  the  continuation  of  rainfall  until  all 
natural  storage  is  satisfied; 

and, 
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5.  the  period  between  termination  of  rainfall 
and  the  first  phase. 

The  five  phases  provide  a  framework  for  analyzing 
the  influence  that  soil  conservation  practices  have  had 
on  surface  runoff  and  baseflow. 

The  first  phase  determines  soil  moisture  just  prior 
to  rainfall  and  influences  baseflow.  A  high  soil 
moisture  content  before  rainfall  produces  increased 
surface  runoff.  The  influence  of  initial  soil  moisture 
on  surface  runoff  ceases  when  the  soil  becomes  saturated 
during  rainfall  (Hawley  et  al.,  1983).  Baseflow  depends 
on  available  moisture  within  the  watershed  and 
groundwater.  Prolonged  periods  between  rainstorms 
decrease  available  soil  moisture  within  the  watershed 
and  deplete  groundwater  reserves  that  in  turn  decrease 
magnitudes  of  baseflow  prior  to  rainfall.  Soil 
conservation  practices  increase  available  moisture 
within  a  watershed  by  retaining  rainfall  within  the 
watershed  and  by  releasing  water  into  drainage  channels 
at  a  reduced  rate. 

Hoyt's  second,  third,  and  fourth  phases  represent 
gradual  accumulation  of  precipitation  and  the  filling  of 
storage  reservoirs.  Strip  cropping  and  contour  farming 
influence  both  baseflow  and  surface  runoff  by  decreasing 
exposed  soil  and  effective  slope  length,  increasing 
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storage  capacity,  and  improving  soil  infiltration 
capacity.  Overall,  soil  conservation  practices  reduce 
surface  runoff  during  this  phase. 

During  Hoyt's  fifth  phase  baseflow  continues  while 
surface  runoff  ends.  A  temporary  new  higher  baseflow 
may  occur  if  infiltration  has  been  sufficient  to  reach 
the  zone  of  saturation.  Soil  conservation  practices 
increase  the  infiltration  rate  enhancing  the  possibility 
of  increased  baseflow.  The  fifth  phase  may  recycle  into 
the  first  phase. 

Artificial  drainage  by  drain  tiles  influences  the 
runoff  process  by  reducing  groundwater  levels  and 
increasing  the  rate  of  subsurface  flow  into  drainage 
channels.  Agricultural  watersheds  utilizing  drain  tiles 
and  other  forms  of  artificial  drainage  tend  to  have 
gentle  slopes  with  minimal  relief;  under  those 
conditions  saturated  soils  hinder  cultivation.  The 
relief  and  stream  density  of  the  Galena  watershed  is 
sufficiently  great  that  little  artificial  drainage  is 
practiced.  Records  pertaining  to  artificial  drainage 
within  the  Galena  watershed  do  not  exist.  Conversations 
with  U.S.G.S.  personnel  and  aerial  photograph  analysis 
during  this  study  period  indicate  no  artificial  drainage 
within  the  watershed. 


There  is  an  extensive  body  of  literature  describing 
the  physical  geography  of  southwestern  Wisconsin  during 
the  first  half  of  the  20th  century.  Blanchard  (1924), 
Martin  (1932),  Trewartha  (1940),  and  Knox  (1972;  1977  ; 
and  1987)  provide  the  interested  reader  detailed 
discussions  about  the  vegetation  and  surface  hydrology 
of  the  Driftless  area  since  the  early  1800s.  The  intent 
of  this  study  is  not  to  conduct  an  extensive  literature 
review  of  previous  studies,  but  to  focus  on  a 
quantitative  runoff  analysis.  The  following  is  a  brief 
discussion  of  key  literature  highlighting  aspects  most 
relevant  to  this  study. 

Blanchard  (1924)  and  Trewartha  (1940)  described  the 
vegetation  encountered  by  the  first  settlers  in 
southwestern  Wisconsin.  The  topography  and  climate  of 
the  Driftless  area  favored  forests  consisting  mainly  of 
oak  that  "fringed  the  streams  and  occupied  the  valleys" 
with  upland  prairies  consisting  of  "thick,  tough,  sod" 
(Blanchard,  1924,  pp.  68  -  69).  The  dense,  thick 
vegetation  prior  to  the  introduction  of  agriculture 
provided  a  very  stable  land  cover  with  good  infiltration 
minimizing  surface  runoff. 

A  dominant  crop  of  the  southern  Driftless  area  from 
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the  earliest  settlers  through  present  times  is  corn. 

The  influx  of  corn  and  other  cultivation  radically 
changed  the  surface  hydrology  by  reducing  the  natural 
vegetation,  decreasing  infiltration,  and  greatly 
increasing  the  portion  of  total  runoff  that  was  surface 
runoff. 

Many  previous  studies  concerning  changing  surface 
hydrology  in  the  Midwest  focus  on  historical 
sedimentation  to  illustrate  changes  in  land  use  since 
the  early  1800s.  Knox  (1972;  1977;  1987)  cites  the 
destruction  of  natural  vegetation  during  cultivation  of 
uplands  and  valley  sides  in  southwestern  Wisconsin  since 
1830  as  apparent  causes  for  increased  magnitudes  and 
frequency  of  peak  discharges,  enlarged  channel  cross 
sections  in  headwater  and  tributary  reaches,  and 
considerable  vertical  accretion  of  sediment  on 
floodplains  in  downstream  valleys. 

Conversion  of  natural  land  to  agricultural  land  use 
since  the  1830s  caused  a  three-to-f ive-fold  increase  in 
the  magnitudes  of  floods  in  the  Platte  watershed  of 
southwestern  Wisconsin,  a  watershed  adjacent  to  the 
Galena  watershed,  indicating  a  significant  increase  in 
surface  runoff  (Knox,  1977)  .  Estimates  of  runoff 
changes  between  the  1830s  and  the  1960s  in  the  Bear 
Branch  tributary  of  the  Platte  River  in  southwestern 
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Wisconsin  indicate  a  tripling  of  peak  flows  from  surface 
runoff  of  moderate  magnitude  storms.  Prior  to  the 
1950s,  increases  in  peak  flows  probably  exceeded  the 
three-fold  magnitude  of  the  1960s  due  to  the  lack  of 
agricultural  conservation  practices  such  as  strip 
cropping  and  contour  plowing  (Knox,  1977).  Knox  (1977) 
identified  a  sequence  of  hydrologic  responses  from  the 
conversion  of  natural  land  cover  to  agricultural  land 
use  and  cultivation  practices:  1830s  -  1860s,  low 
magnitude  disturbance;  1870s  -  1940s,  maximum 
disturbance;  1950s  -  1970s,  moderate  disturbance.  The 
period  from  1950s  -  1970s  extends  into  the  1980s. 
Hydrologic  disturbance  appears  to  subside  with  time  in 
this  later  period  due  to  the  implementation  of 
agricultural  conservation  practices. 

Trimble  and  Lund  (1982)  compared  agricultural 
practices  of  the  1930s  and  the  1970s  to  evaluate  the 
influence  of  agricultural  conservation  practices  on  soil 
erosion  in  the  Coon  Creek  Basin,  Wisconsin.  The  area 
surrounding  Coon  Creek  in  1934  consisted  of  rectangular 
fields  on  moderate  to  steep  slopes,  poor  crop  rotations, 
removal  of  crop  residues,  insufficient  manuring, 
nutrient  depletion,  lack  of  cover  crops,  and  very  active 
erosion.  The  same  area  in  the  mid  1970s  portrayed 
contour  plowing,  strip  cropping,  long  crop  rotations, 
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incorporation  of  crop  residues  into  the  soil,  manuring, 
and  cover  crops.  The  improvement  in  agricultural 
practices  within  the  Coon  Creek  Basin  has  increased  soil 
infiltration  capacities  and  converted  what  would  have 
been  a  significant  portion  of  surface  runoff  into 
subsurface  flow,  thereby  reducing  soil  erosion.  Trimble 
and  Lund  discount  a  climatic  influence  believing  that 
climatic  cycles  superimpose  perturbations  on  overall 
trends  caused  by  agricultural  conservation  practices. 

Sartz  (1976)  monitored  runoff  response  to  the 
influence  of  different  land  uses  at  the  Coulee 
Experimental  Forest  near  LaCrosse,  Wisconsin.  Sartz 
found  that  runoff  from  cultivated  land  was  ten  times 
greater  than  a  lightly  grazed  pasture  and  approximately 
20  times  greater  than  runoff  from  a  forested  pasture. 
Agricultural  conservation  practices  were  not  evaluated. 

Sartz  (1974)  studied  the  influence  that 
non-cultivation  agricultural  practices  have  on  runoff  by 
comparing  grazed  watersheds  in  southwestern  Wisconsin. 
Runoff  magnitudes  decreased  several  times  within  only 
two  years  after  grazing  was  discontinued  in  one 
watershed,  illustrating  the  wide  variety  of  agricultural 
conservation  practices  available  to  decrease  runoff. 

C.M.  Adamson  (1974)  studied  two  similar  small 
watersheds  in  Australia  focusing  on  the  influence  that 
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soil  conservation  practices  have  on  surface  runoff. 

Both  watersheds  consisted  of  undisturbed  natural  prairie 
at  the  beginning  of  the  study.  Adamson  maintained  one 
watershed  at  its  original  condition  while  contour 
plowing  and  fertilization  were  applied  to  the  other 
watershed.  Results  indicate  a  reduction  in  surface 
runoff  of  74%  in  the  treated  watersned  during  the  21 
year  study.  The  reduction  was  attributed  to  increased 
infiltration  and  water  storage  capacity  from  contour 
plowing. 

SUMMARY 

The  runoff  process  consists  of  many  interrelated 
factors  influenced  by  land  use.  Soil  conservation 
practices  influence  runoff  from  agricultural  land  by 
reducing  surface  runoff  and  increasing  baseflow.  Poor 
documentation  of  historical  events  complicate  an 
accurate  description  of  surface  runoff  response  to 
changing  land  use. 
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CHAPTER  3:  DESCRIPTION  OF  THE 
GALENA  WATERSHED 

The  Galena  watershed  is  within  the  Driftless  area 
of  southwestern  Wisconsin  and  northwestern  Illinois. 

The  study  area  is  upstream  of  the  Buncombe,  Wisconsin 
U.S.  Geological  Survey  gaging  station  and  lies  primarily 
in  Lafayette  County,  Wisconsin  and  a  small  area  within 
Jo  Daviess  County,  Illinois  (Figure  1).  The  U.S.G.S. 
operates  a  continuous  recording  discharge  gage  located 
at  the  SW  1/4,  Section  33,  TIN,  R1E,  4th  P.M. ,  in 
Lafayette  County,  Wisconsin.  Operation  of  the 
continuous  recording  discharge  gage  began  in  December, 
1939.  The  medium-size  Galena  watershed  is  323  square 
kilometers  in  area.  Terrain  within  the  watershed  is 
highly  dissected  with  local  relief  ranging  from  5-15 
meters  per  square  kilometer  (m/km2)  in  the  northern 
portion  to  70  -  100  m/km2  the  steep  southern  valleys. 
Urbanization  within  the  watershed  is  virtually 
non-existent  with  only  a  few  small  towns  and  numerous 
clusters  of  houses.  Initial  land  use  in  the  watershed 
was  mining  until  the  1830's  when  agriculture  became  the 
predominant  land  use  (Knox,  1987).  The  Galena  watershed 
was  chosen  for  study  because  of  the  availability  of  good 
historical  records,  aerial  photographs,  and  because 
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significant  change  in  agricultural  land  use  has  occurred 
within  the  last  fifty  years. 

The  surficial  bedrock  of  this  area  is  Ordovician 
age  and  consists  mainly  of  Galena  dolomite  with  local 
exposures  of  Maquoketa  shale  (Agnew,  1963;  Klemic  and 
West,  1964;  Mullens,  1964).  The  drainage  pattern  is 
dendritic  and  rectangular  showing  strong  geologic 
influence  through  joint  controlled  valleys  and  lineated 
channel  patterns.  The  shape  of  this  watershed  is 
elongated  with  the  primary  axis  oriented  slightly 
northwest  to  southeast.  The  orientation  decreases  the 
effectiveness  of  storm  movement  to  magnify  floods 
because  most  storms  track  from  the  south  and  west 
to  the  north  and  east  (Bryson,  1966) . 

There  are  three  soil  associations  within  the  Galena 
watershed  associated  primarily  with  loess  deposited 
during  the  late  Wisconsin  (Watson,  1966,  Knox,  1987) . 

The  Tama-Ashdale  Association  represents  dark-colored, 
deep  soils  of  the  limestone  uplands  found  mostly  on  the 
numerous  broad  ridge  tops  and  adjoining  slopes.  Tama 
soils  are  approximately  125  centimeters  thick  while  the 
Ashdale  soils  average  75  centimeters.  Both  Tama  and 
Ashdale  soils  belong  to  the  Mollisol  soil  order.  Tama 
soils  have  an  infiltration  rate  between  2.0  -  6.4 
centimeters  per  hour  (cm/hr)  throughout  the  soil.  The 
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Ashdale  soils  have  an  infiltration  rate  of  2.0  -  6.4 
cm/hr  from  the  land  surface  to  30  centimeters  depth,  0.5 
-  2.0  cm/hr  from  30  -  100  centimeters  depth  and  0.13  - 
0.5  cm/hr  from  100  to  approximately  125  centimeters 
depth  (until  the  limestone  parent  material) .  The 
Tama-Ashdale  association  covers  approximately  40%  of  the 
watershed. 

The  Fayette-Palsgrove  association  represents 
light-colored,  deep  soils  also  found  on  the  limestone 
uplands.  The  association  also  represents  approximately 
40%  of  the  watershed.  Small  piles  of  gravelly  and  stony 
waste  material  from  lead  and  zinc  mines  ranging  in  size 
from  2-20  acres  are  found  on  the  the  Fayette-Palsgrove 
Association.  Both  Fayette  and  Palsgrove  soils  are  75  - 
125  centimeters  thick  and  belong  to  the  Alfisol  soil 
order.  The  Fayette  soils  have  an  infiltration  rate  of 
2.0  -  6.4  cm/hr  throughout  the  soil  profile.  The 
Palsgrove  soils  have  the  same  infiltration  rate  as  the 
Ashdale  soils. 

The  Dubuque-Sogn  Association  represents 
light-colored  moderately  deep  and  dark  colored  shallow 
soils.  Dubuque  soils  are  38  -  75  centimeters  thick 
while  the  Sogn  soils  are  the  thinnest  soils  in  the 
watershed,  averaging  30  centimeters  or  less.  Dubuque 
soils  belong  to  the  Alfisol  soil  order  and  have  an 
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I  infiltration  rate  of  2.0  -  6.4  cm/hr  from  the  surface  to 

I  30  centimeters  depth,  0.5  -  2.0  cm/hr  from  30  -  75 

centimeters  depth,  and  0.13  -  0.5  cm/hr  from  75 
|  centimeters  depth  to  the  limestone  parent  material. 

Sogn  soils  are  the  least  fertile  agricultural  soils  and 
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they  are  found  mostly  on  steep  side  slopes  bordering 
stream  valleys.  They  are  very  droughty,  lack  space  for 
root  development,  and  they  are  stony  in  many  places. 

Sogn  soils  belong  to  the  Entisol  soil  order  and  they 
have  an  infiltration  rate  of  2.0  -  6.4  cm/hr  throughout 
their  profiles.  The  Dubuque-Sogn  Association  covers 
approximately  20%  of  the  watershed. 


This  study  utilizes  three  data  sets: 
precipitation,  stream  runoff,  and  land  use. 

Precipitation  data  consist  of  records  from  four  weather 
stations  surrounding  the  Galena  watershed.  Hydrographs 
of  the  Galena  River  provide  necessary  information  to 
estimate  surface  runoff  and  baseflow  data.  Analyses  of 
aerial  photographs  provide  land  use  data.  All  data  sets 
span  the  study  period  with  minimal  missing  data. 
Statistical  analyses  were  conducted  to  evaluate  linkages 
between  precipitation,  runoff,  and  land  use  data  sets. 

PRECIPITATION 

Precipitation  occurring  at  weather  stations  located 
at  Cuba  City,  Darlington,  and  Platteville,  Wisconsin, 
and  Galena,  Illinois  (Figure  2)  was  areally  weighted  by 
the  Thiessen  Polygon  method  (Dunne  and  Leopold,  1978)  to 
estimate  rainfall  occurring  within  the  watershed.  The 
Cuba  City  station  is  the  most  important  using  the  Thies¬ 
sen  method  and  is  the  only  hourly  recording  station  of 
the  four.  This  study  only  considered  rainfall  occurring 
from  April  through  October  to  avoid  snowmelt  and  cold 


weather  influences.  To  ensure  the  occurrence  of  stream 
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Figure  2.  The  spatial  location  of  the  four  weather 
stations  used  In  this  study  to  deteraine  the 
precipitation  and  temperature  data. 
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runoff  following  precipitation,  a  decision  was  made  to 
include  in  the  analyses  only  those  events  for  which  at 
least  2.50  centimeters  of  precipitation  occurred  within 
24  hours.  A  total  of  145  storms  equalling  2.50  cm  have 
occurred  within  the  study  period  1940-1987.  After 
adjusting  to  ensure  homogeneity  of  the  climate  data  base 
(further  discussion  in  Chapter  Five),  83  storms  were 
selected  for  analyses.  The  selected  storms  were 
subdivided  into  three  time  periods  based  on  the  changing 
percentage  of  soil  conservation  practices  obtained  from 
land  use  analysis.  The  three  resulting  time  periods 
are:  1940  -  1950,  initial  implementation  of  soil 

conservation  practices;  1952  -  1968,  moderate  increase 
in  soil  conservation  practices  observed;  and  1969  - 
1987,  representing  a  maximum  level  for  application  of 
soil  conservation  practices.  Statistical  analyses 
showed  that  storm  intensities  were  similar  between  the 
three  time  periods. 

SURFACE  RUNOFF/BASEFLOW 

The  estimated  percentages  of  surface  runoff  and 
baseflow  for  each  storm  result  from  quantifying 
components  of  corresponding  storm  hydrographs. 
Hydrographs  of  the  Galena  River  taken  at  Buncombe, 
Wisconsin,  for  the  selected  storms  were  obtained  from 
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the  U.S.G.S.  Water  Resources  Division.  Surface  runoff 
duration  for  each  storm  resulted  from  a  visual 
inspection  of  each  hydrograph,  denoting  the  interception 
of  the  rising  limb,  the  peak  stage  position  on  the 
hydrograph,  and  the  inception  of  the  asymtotic  position 
of  the  hydrograph.  The  surface  runoff  volume  is 
estimated  by  measurement  of  the  hydrograph  area  above  a 
straight  line  connecting  the  inception  points  on  the 
rising  limb  and  the  asymtotic  recession  limb.  Base  flow 
is  the  area  of  the  hydrograph  below  that  line.  This 
method  is  not  exact,  but  it  is  consistent  throughout 
this  analysis  and  therefore  it  is  reasonably  objective. 

Each  data  base  describing  surface  runoff  and 
baseflow  consists  of  a  series  of  (X,Y)  coordinates.  The 
X  axis  represents  time  in  hours  and  the  Y  axis 
represents  discharge.  A  volume  representing  surface 
runoff  and  baseflow  for  each  storm  results  from 
processing  storm  data  bases  through  a  Fortran  program 
(Appendix  1) . 

LAND  USE 


This  study  uses  aerial  photographs  to  analyze 
changing  land  use  during  the  study  period.  Aerial 
photographs  taken  at  irregular  intervals  by  the 
Agriculture  Stabilization  and  Conservation  Service 
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(A.S.C.S.)  and  the  U.S.G.S.  began  in  1937.  Photographs 
used  in  this  study  are  from  1937,  1951,  1955,  1962, 

1968,  1976,  and  1985. 

Percentages  of  land  areas  in  various 
classifications  were  estimated  by  observing  land  uses 
along  twelve  east/west  transects  arranged  from  north  to 
south  within  the  Galena  watershed  (Figure  3) .  Distance 
between  transects  is  approximately  two  kilometers  to 
ensure  an  accurate  account  of  land  use  in  small  order 
tributaries  and  headwaters.  Transect  locations  were 
determined  using  permanent  land  marks  such  as  road 
junctions  and  bridges.  An  opsometer  was  used  to  measure 
the  four  classifications  of  land  use.  Land  use 
percentages  resulted  from  the  total  length  of  the 
transects. 

Photograph  clarity  restricts  land  use 
classification  in  this  study  to  cultivated, 
non-cultivated,  and  forest.  Analyses  of  the  culti¬ 
vated  land  determined  the  percentage  of  land  associated 
with  soil  conservation  practices,  primarily  strip 
cropping  and  contour  farming. 

STATISTICAL  ANALYSIS 

A  comparison  of  rainfall  intensities  between  the 
three  time  periods  was  conducted  by  using  analysis  of 
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Flgure  3.  The  epatlai  Location  -*  the  twelve  tranaecte 
used  to  evaluate  land  use. 
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variance  and  other  descriptive  statistics.  These 
analyses  showed  storm  intensities  between  the  three 
periods  were  not  significantly  different  for  the 
adjusted  data. 

Multiple  regression  analyses  were  applied  to 
determine  relationships  and/or  contributions  of  eight 
hydrologic  variables  that  influence  stream  runoff. 
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CHAPTER  5:  DATA  ANALYSIS 


INTRODUCTION 


The  objectives  of  this  chapter  are: 

1.  to  explain  the  criteria  used  in  selecting 
analyzed  storms; 

2.  to  analyze  land  use  within  the  Galena 
watershed  during  the  study  period; 

3.  to  analyze  storm  and  runoff  characteristics 
that  could  influence  the  percentage  of 
surface  runoff  occurring  during  this  study 
period; 

and, 

4.  to  compare  storm  and  runoff  characteristics 
between  three  periods  of  soil  conservation 
intensity. 

This  chapter  analyzes  and  compares  characteristics 
of  runoff  events  determined  from  hydrographs  and  storms 
from  1940  through  1987  to  evaluate  hydrologic  response 
to  increased  utilization  of  soil  conservation 
practices.  Aerial  photograph  analysis  defined  three 
periods  of  varying  degrees  of  soil  conservation  activity 
within  the  Galena  watershed  during  the  study  period. 
Aerial  photographs  were  used  to  classify  land  use  within 
the  Galena  watershed  during  this  study  period  into 
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percent  cultivated,  non-cultivated,  and  forested. 

Further  analysis  determined  the  percentage  of  cultivated 
land  utilizing  soil  conservation  practices.  Multiple 
regression  analysis  was  then  used  to  evaluate 
relationships  between  storm  characteristics  and  runoff 
hydrographs  for  the  study  period. 

STORM  SELECTION 

Storm  rainfalls  were  selected  from  the  records  of 
four  weather  stations  within  or  surrounding  the  Galena 
watershed.  The  stations  are  located  at  Meeker's 
Grove/ Cuba  City,  Darlington,  and  Platteville,  Wisconsin, 
and  Galena,  Illinois  (Figure  2) .  The  only  hourly 
recording  station  among  the  four  is  the  Meeker's 
Grove/Cuba  City  station  located  on  the  the  western  edge 
of  the  watershed.  That  record  was  used  to  estimate 
storm  durations  of  the  other  three  stations.  The  bias 
introduced  by  estimating  storm  durations  for  the  other 
three  stations  is  thought  to  be  insignificant  because 
the  Meeker's  Grove/ Cuba  City  station  is  strongly 
dominant  (87%/74%  from  the  Thiessen  Polygon  method)  in 
the  allocation  of  rainfall  on  the  watershed.  Station 
histories  from  1940  to  1988  indicate  that  all  stations 
moved  slightly  on  various  occasions,  but  only  on> 
significant  move  occurred  when  the  station  located  at 
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Meeker's  Grove  moved  approximately  6.5  kilometers 
southwest  to  Cuba  City  in  June  1951.  It  is  not  possible 
to  calibrate  the  effects  of  the  moves  for  any  of  the 
stations  due  to  the  lack  of  overlap  in  monitoring  at  the 
stations.  The  records  of  the  four  weather  stations  are 
very  good  with  minimal  missing  data.  Values  for  the 
occasional  missing  data  were  estimated  using  regression 
equations  related  to  data  from  the  other  stations  (Dunne 
and  Leopold,  1978) . 

Because  available  weather  stations  are  few  in 
number  and  unevenly  distributed  in  relation  to  the 
Galena  study  watershed,  the  Thiessen  Polygon  method 
(Dunne  and  Leopold,  1978)  was  applied  to  determine  the 
relative  watershed  allocations  of  rainfall  from  the  four 
stations.  Thiessen  polygon  allocations  show  the 
respective  influences  of  the  precipitation  events  at  the 
four  stations  on  watershed  runoff  from  1940  to  June  1951 
are:  Meeker's  Grove  -  87%,  Darlington  -  4%,  Platteville 
-  4%,  and  Galena  -  5%,  and  from  June  1951  to  1988:  Cuba 
City  -  74%,  Darlington  -  10%,  Platteville  -  14%,  and 
Galena  -  2%.  The  Thiessen  Polygon  allocations  were 
analyzed  by  correlating  the  rainfall  of  the  four 
stations  to  the  gage  heights  recorded  at  the  Buncombe, 
WI,  discharge  station.  This  analysis  linked  a  physical 
property,  the  discharge  gage  height,  to  the  spatial 
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distribution  of  the  four  stations.  The  resulting 
correlations  supported  the  Thiessen  allocations. 

Assessment  end  comparison  of  the  influence  of 
changing  land  use  on  the  relative  fractions  of  surface 
versus  baseflow  runoff  for  the  three  time  periods 
required  selecting  storms  that  were  as  similar  as 
possible  in  amount  and  duration  of  precipitation.  The 
expected  two-year  probability  rainfall  for  the  study 
watershed  ranges  from  approximately  2.5  cm  in  30  minutes 
to  approximately  7.5  cm  in  24  hours  (Dunne  and  Leopold, 
1978,  pp.  58  -  63).  Based  on  predicted  intensities,  the 
need  to  analyze  storms  that  occur  frequently,  and  to 
ensure  a  large  sample  size  of  storms  that  have 
sufficient  intensity  to  produce  significant  surface 
runoff,  a  minimum  storm  intensity  of  2.5  centimeters 
occurring  within  24  hours  was  defined  as  a  minimum 
threshold  for  inclusion.  Only  storms  occurring  from 
April  through  October  were  considered  to  avoid  the 
influence  of  snowmelt  and  frozen  ground.  The  threshold 
definition  resulted  in  identification  of  145  storms  from 
1940  to  1987. 

The  Thiessen  method  was  applied  with  the  other 
three  stations  on  the  initial  sample  of  145  storms. 
Analysis  of  variance  was  then  conducted  on  the  initial 
145  storm  sample  precipitations  and  durations  among  the 
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three  time  periods  (Table  1) .  1  All  analyses  indicated 
similarity  among  time  period  precipitations  magnitudes. 
However,  a  significant  difference  (0.03  probability 
level)  occurred  between  the  1940s  storm  duration  (mean 
of  19.12  hours)  and  the  other  two  time  period  storm 
durations  (means  of  14.00  and  13.90  hours).  The 
influence  of  the  slightly  longer  storm  duration  for  the 
1940s  is  believed  to  be  minimal.  Overall,  statistical 
analyses  strongly  indicate  that  climatic  influences  are 
similar  throughout  the  study  period  and  should  not 
influence  runoff  analyses.  However,  to  ensure 
homogeneity  of  the  climatic  influences  the  initial 
sample  of  145  storms  was  adjusted  resulting  in  a  final 
sample  set  of  83  storms. 

The  initial  sample  of  145  storms  were  plotted  on  an 
X/Y  coordinate  system  with  the  X  axis  representing 
duration  and  the  Y  axis  representing  the  amount  of 
precipitation  (Figure  4) .  Extreme  outliers  unique  to 
only  one  of  the  three  groups  were  eliminated  to  promote 
similarity  among  the  groups.  Also  eliminated  were 
storms  occurring  within  five  days  of  a  previous  storm 
because  they  created  a  circumstance  whereby  the 
antecedent  precipitation  might  greatly  bias  surface 
runoff  (Mockus,  1972,  and  phone  conservation  with  John 
Milligan,  Soil  Conservation  Service,  Madison,  WI,  June 
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NEAR  STANDARD  35*  C.I.  F  -  RATIO  PROBABILITY 

DEVIATION  FOR  PERIOD  FOR  THREE  PERIODS 


STORM  PRECIPITATION 
TIME  PERIOD  I 
i  54(/  -  !  95  i 

(IN  CM) : 

4.76 

1  Cl 

I  .71 

0  OA  _  C  On 

0.8! 

0.45 

TIME  PERIOD  2 
]  352  -  !3E.§ 

C  Or. 

v* .  iV 

I 

1  qf.  _  c  Of. 

TIME  PERIOD  3 

1863  -  1387 

5.18 

-*>  1  1 
L.  1  1 

4.80  -  5.90 

STORM  DURATION  (IN  1 
TINE  PERIOD  1 

•.Gift  _  1GC1 

!?•»•/  1  JJ  1 

TOURS) : 

1  0  !*) 

!  J  .  !  i 

11.64 

17.1  -  21.4 

3.70 

0.03 

TINE  PERIOD  2 

1352  -  1368 

14.01 

7.81 

13.2  -  17.3 

TINE  PERIOD  3 

1363  -  1387 

13.31 

\  1  40 

1  1  ,  H'J 

13.1  -  17.3 
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1989) .  The  resulting  83  storms  for  the  three  time 
periods  include  27  storms  occurring  from  1940  through 
1950,  29  storms  occurring  from  1952  through  1968,  and  27 
storms  occurring  from  1969  through  1987  (Figure  5) . 
Analysis  of  variance,  with  a  research  hypothesis  that 
u1=u2=u3  for  storm  durations  and  precipitation 
magnitudes  indicate  similarity  between  the  three  groups 
(Table  2) .  Furthermore,  a  plot  of  storm  precipitation 
verses  time  shows  no  statistically  significant  trend  and 
the  resulting  regression  equation  is  Y  =  0.007 IX  -  9.36 
(Figure  6  and  Table  3) .  An  additional  indication  of 
similarity  results  from  comparing  precipitation  to  storm 
duration  for  each  storm.  The  resulting  ratio  gives  a 
relative  indication  of  storm  intensity.  Average  ratios 
and  statistical  analyses  for  the  three  time  groups 
indicate  similarity  (see  Tables  3  and  4) . 

It  is  virtually  impossible  to  precisely  select 
storms  that  are  similar  in  antecedent  conditions  due  to 
the  wide  variety  and  spatial  variations  for  antecedent 
conditions  within  a  watershed.  With  a  large  number  of 
storms  in  each  group  the  wide  variance  of  antecedent 
conditions  is  represented  in  each  group  increasing  the 
similarity  among  the  three  groups  and  decreasing  the 
bias  from  antecedent  conditions.  Antecedent 
precipitation  is  analyzed  with  other  independent 
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Figure  5.  A  plot  of  the  83  storm  intensities  after 
eliminating  outliers  and  those  storms  not  meeting 
antecedent  precipitation  criteria. 
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TABLE  2  STATISTICS  OH  THE  83  STORK  SAMPLE 


MEAN  STANDARD  95*  C.I.  F  -  RATIO  PROBABILITY 

DEVIATION  FOR  PERIOD  FOR  THREE  PERIODS 

STORK  PRECIPITATION  (IN  CM):  0.8?  0.42 

TIKE  PERIOD  1 


1940  -  1951 

4.43 

i  .08 

0  00  -  4  0*7 

•J  .JJ  n.’Ji 

Ttur  ocomr.  7 

tine  <  4 

1952  -  1958 

4  83 

1 

1  .  L’J 

1  J1  -  c  oc 

H.*t(  -j  .  L-J 

TIKE  PERIOD  3 

1969  -  1987 

4.56 

\  1  1 
t  .  (  1 

4.13  -  5.05 

STORK  DURATION  (IN 
TIKE  PERIOD  I 

HOURS); 

0.44 

0.64 

1940  -  1951 

15.74 

9.57 

12.48  -  19.00 

TIKE  PERIOD  2 

1952  -  1968 

13.52 

7.51 

10.25  -  16.75 

TIKE  PERIOD  3 

1969  -  1987 

14.48 

9.37 

11.00  -  17.75 

RATIO  (PRECIPITATION /DURAT ION): 

0.52 

0.60 

TIKE  PERIOD  1 

1940  -  1951 

0.42 

A  00 

V  .  L'J 

0.29  -  0.54 

TIKE  PERIOD  2 

1952  -  1968 

0.48 

0.31 

0.35  -  0.60 

TIKE  PERIOD  3 

1969  -  '987 

0.51 

0  41 

0.38  -  0.64 
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TABLE  3  STATISTICS  ON  STORM  6R0UPS 

A:  STORM  DISTRIBUTION  BY  PRECIPITATION  AMOUNTS  (IN  CM) 


PRECIPITATION 

AMOUNT: 

2.50  -  4  00 

4 . 00 

-  6.00 

6.00  - 

8.00 

TIME  PERIOD  I 
1340  -  1351 

\  * 
i  i. 

12 

3 

1 

TIME  PERIOD  2 
1352  -  1963 

1 

t 

17 

5 

TIME  PERIOD  3 
1963  -  1987 

9 

14 

4 

B:  STORM  DISTRIBUTION  BY  MONTHS 

MONTHS: 

APRIL  MAY 

■JUNE 

JULY 

AUGUST 

SEPTEMBER  OCTOBER 

TIME  PERIOD  1 
1940  -  1351 

3  4 

8 

1 

6 

4  1 

TIME  PERIOD  2 
1352  -  1969 

2  3 

5 

6 

6 

4  3 

TIME  PERIOD  3 
1363  -  1987 

4  2 

2 

5 

0 

J 

4  1 

C:  REGRESSION  ANALYSIS  -  AVERAGE  MONTHLY  PRECIPITATION  (APR  -  OCT)  VS  TIME: 

CORRELATION  OF  THE  AVERAGE  MONTHLY  PRECIPITATION  VS  TIME;  0.059 

REGRESSION  EQUATION:  Y  =  0.0071X  -  9.36 

STANDARD  t  -  VALUE  PROBABILITY  F  -  RATIO  PROBABILITY 

ERROR  (X  =  0)  LEVEL  LEVEL  R  SQUARED 

0.0! 9  0.40  0.69  0.16  0.69  0.0034 


TABLE  i  STOSS  ANS  HYCROGSAPH  SATA 
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variables  utilizing  regression  analysis. 

LAND  USE 

Land  use  was  evaluated  by  estimating  from  aerial 
photographs  the  percentages  of  a  given  land  use  observed 
along  line  transects  drawn  across  the  Galena  watershed 
at  specified  intervals.  Twelve  transects  spaced 
approximately  two  kilometers  apart  and  oriented  in  a 
predominantly  east/west  direction  provided  adequate 
spatial  coverage  to  ensure  a  representative  sample  of 
lower  order  streams  and  headwaters  (Figure  3) .  Scale 
differences  between  photograph  series  were  not  a  bias  as 
this  analysis  is  based  on  relative  percentages.  The 
large  scale  of  the  aurial  photographs  allowed  three  land 
use  classifications:  cultivated,  non-cult ivated,  and 
forested.  A  further  division  of  the  cultivated  land 
into  that  with  soil  conservation  practices  verses  that 
with  no  soil  conservation  practices  was  also 
undertaken.  Classification  of  row  crops  and  most  other 
small  scale  land  uses  were  not  possible  due  to  the  scale 
and  clarity  of  the  aerial  photographs.  Existing 
agricultural  statistics  describe  agricultural  land  use 
on  a  county  level  prohibiting  a  specific  analysis  of  the 
Galena  watershed.  All  known  aerial  photography  taken 
during  this  study  period  was  analyzed.  These  included 
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photographs  for  the  following  years:  1937,  1951,  1955, 
1962,  1968,  1976,  and  1985. 

Table  5  summarizes  the  changing  land  use  within  the 
Galena  watershed  from  1937  through  1985.  All 
photographic  series  except  the  1951  and  1968  series  were 
taken  from  late  August  to  late  October  (see  Table  5) . 

The  possible  bias  introduced  from  analyzing  photographs 
from  different  seasons  is  minimal  due  to  the  broad  land 
use  classification  categories  and  to  the  ease  in 
identifying  soil  conservation  practices  such  as  strip 
cropping  and  contour  plowing.  One  possible  bias  from 
analyzing  photographs  from  different  seasons  might  occur 
as  leaf  cover  appears  greater  in  autumn  than  spring. 

The  fluctuating  percentage  of  forested  land  may  be 
biased,  although  the  1985  series  does  not  conform  to 
this  bias.  The  small  percentage  of  forested  land 
throughout  the  study  period  decreased  the  significance 
of  this  possible  bias. 

Land  use  varies  from  north  to  south  within  the 
Galena  watershed.  The  northern  portion  of  the  Galena 
watershed  has  less  relief  and  is  characterized  by  large, 
gently  sloping,  cultivated  fields  utilizing  less  soil 
conservation  practices  than  in  fields  of  the  steeper 
southern  portion.  The  majority  of  cultivation  occurs  in 
the  northern  portion  while  the  southern  portion  has  the 


50 


majority  of  forested  land. 

The  1937  and  1951  aerial  photographs  indicate 
minimal  use  of  soil  conservation  practices  (Table  5  and 
Figure  7) .  Cultivated  fields,  especially  those  in  the 
southern  part  of  the  watershed,  indicate  extensive 
gullying  and  soil  erosion.  Forests  appear  sparse 
compared  to  that  shown  on  later  aerial  photograph 
series.  Overall,  the  land  use  from  1937  to  1951 
indicates  poor  conservation  practices. 

From  1955  through  1968  adoption  of  soil 
conservation  practices  increased.  Gullying  appears  less 
extensive  and  headward  erosion  of  channels  appears  to 
subside.  Forests  appear  more  dense.  The  influence  of 
improved  agricultural  technology  is  suggested  during 
this  time  period  as  the  1937  photographs  indicated  corn 
shocks  geometrically  spaced  in  denuded  and  harvested 
fields  while  the  1955  and  later  photographs  are  devoid 
of  any  indication  of  corn  shocks. 

The  latest  period  indicates  a  substantial  increase 
in  the  use  of  soil  conservation  practices.  Gullying 
exists,  but  is  at  its  lowest  activity  level  since  the 
late  1930s  and  1940s  and  is  characterized  by  minimal 
headward  erosion.  The  quality  of  soil  conservation 
practices  appears  to  improve  during  this  period  as 
contour  farming  is  more  precise,  strip  cropping  is  more 
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frequent,  and  the  use  of  buffer  strips  around  fields  is 
more  apparent. 

Soil  conservation  practices  in  the  Galena  watershed 
consist  primarily  of  strip  cropping  and  contour 
farming.  From  1937  through  1962  strip  cropping  and 
contour  farming  were  the  only  conservation  practices 
observed.  Beginning  in  1968,  indications  of  terracing 
and  buffer  strips  around  fields  began  to  appear  in 
addition  to  previous  cultivation  conservation  practices. 

Several  federal  and  state  programs  initiated  since 
the  mid-1960s  appear  to  have  had  minimal  influence  on 
the  significant  increase  in  soil  conservation  within  the 
Galena  watershed  since  the  late  1960s.  In  1965  the 
Agriculture  Stabilization  and  Conservation  Service 
(A.S.C.S.)  initiated  the  Soil  Bank  Program  which  gave 
farmers  federal  subsidies  for  taking  cultivated  land  in 
erosion  prone  areas  out  of  cultivation  for  a  period  of 
10  years.  The  program  exists  today  as  the  Conservation 
Reserve  Program  (CRP) .  The  hydrologic  influence  of  the 
Soil  Bank  and  CRP  programs  on  the  Galena  watershed  is 
insignificant  as  less  than  one  percent  of  the  watershed 
area  participated  in  the  programs  during  the  study 
period  (conversation  with  Lafayette  County  A.S.C.S. 
Executive  Director  Leon  Wolfe,  December,  1989) . 

The  Galena  watershed  received  State  funding  during 
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the  early  1970s  to  construct  terraces  in  highly  erodible 
areas.  Statistics  defining  the  amount  of  terracing 
constructed  within  the  Galena  watershed  do  not  exist. 
Terracing  is  difficult  to  differentiate  from  contour 
farming  on  aerial  photographs,  prohibiting  a  precise 
evaluation  of  its  significance  since  the  early  1970s. 
However,  the  amount  of  terracing  within  the  Galena 
watershed  is  estimated  to  be  less  than  five  percent  of 
the  area  utilizing  soil  conservation  practices 
(conversation  with  Lafayette  County  A.S.C.S.  Executive 
Director  Leon  Wolfe,  December  1989) . 

The  most  probable  factor  influencing  the  increase 
in  soil  conservation  since  the  late  1960s  is  the 
farmer's  awareness  and  understanding  of  soil 
conservation  practices.  A  previous  study  by  Johansen 
(1969)  illustrates  the  rapid  rate  farmers  adopted 
contour  strip  cropping  from  1939  (233  users)  to  1967 
(6402  users)  in  southwestern  Wisconsin  as  they  became 
aware  of  the  methods  and  benefits  of  soil  conservation. 

Farmers  have  also  responded  to  the  landmark  1985 
Federal  Food  Security  Act  which  mandates  the  use  of  soil 
conservation  practices  as  a  prerequisite  to  receiving 
Federal  subsidies.  Federal  subsidies  are  extremely 
important  to  most  farmers  and  may  account  for  a 
significant  portion  of  their  financial  income. 
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Lafayette  County  A.S.C.S.  Executive  Director  Leon  Wolfe 
estimates  that  85-90%  of  the  cultivated  land  within  the 
Galena  watershed  will  involve  soil  conservation 
practices  by  1995  as  a  consequence  of  the  Federal  Food 
Security  Act. 

ESTIMATED  SURFACE  RUNOFF  AND  BASEFLOW 

Hydrograph  data  representing  runoff  from  the  83 
sample  storms  were  obtained  from  the  U.S.G.S.  Water 
Resources  Division  archives.  The  U.S.G.S.  gaging 
station  located  on  the  Galena  River  at  Buncombe, 
Wisconsin  operated  from  1940  until  July  1967  with  a 
Friez  FA-3  water-stage  recorder  which  produced  a  paper 
strip  chart  indicating  the  actual  water  stage  of  the 
river  in  feet  (U.S.G.S.  Station  Analysis,  1940  -  1987) . 

A  rating  curve  converts  the  water  stage  to  discharge  and 
is  updated  periodically  as  the  channel  geometry 
changes.  In  July  1967  the  U.S.G.S.  installed  a 
Fischer-Porter  digital  recorder  which  records  the  water 
depth  every  15  minutes  in  a  digital  format  eliminating 
the  paper  strip  chart.  Since  1986,  the  U.S.G.S. 
instituted  a  computerized  system  that  instantaneously 
applies  the  rating  curve  to  the  digital  depth  to  produce 
a  record  of  discharge  in  cubic  feet  per  second  every  15 


minutes. 
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The  U.S.G.S.  prepares  an  annual  station  analysis 
describing  the  water  year  for  each  gage.  The  station 
analysis  defines  periods  when  shifts  in  the  relationship 
between  gage  height  and  discharge  have  occurred.  Shifts 
may  occur  due  to  changes  in  the  channel  geometry  at  the 
gaging  station.  All  applicable  shifts  were  applied  to 
the  hydrographs  analyzed  in  this  study. 

Each  hydrograph  was  guantified  to  obtain  accurate 
percentages  of  surface  runoff  and  baseflow  from  storm 
events.  The  strip  chart  data  or  digital  hydrograph  data 
were  converted  into  discharge  using  the  rating  curves. 
The  discharge  data  were  then  converted  into  an  (X,Y) 
coordinate,  the  X-axis  representing  time  in  hours  and 
the  Y-axis  representing  discharge.  This  process 
produced  a  data  file  containing  (X,Y)  coordinates  that 
were  used  to  quantify  both  surface  runoff  and  baseflow 
for  the  storms  analyzed.  The  portion  of  the  hydrograph 
representing  surface  runoff  was  determined  by  a  visual 
inspection  of  the  hydrograph  data  as  previously 
discussed. 

Data  files  describing  surface  runoff  and  baseflow 
were  processed  through  a  Fortran  computer  program  that 
determines  the  areas  of  irregularly  shaped  polygons 
(Appendix  1) .  The  computer  program  determined  the  areas 
under  the  hydrograph  curve  representing  surface  runoff 
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and  baseflow.  The  area  determined  by  the  Fortran 
program  equated  to  a  volume  as  discharge  (defining  the  Y 
axis  in  CMS)  was  multiplied  by  time  (defining  the  X  axis 
in  hours)  (Table  4) .  The  Fortran  program  converted  the 
hours  to  seconds  and  produced  a  final  output  in  cubic 
meters  which  facilitated  calculating  the  percent  of 
surface  runoff  and  base  flow. 

The  estimated  percentages  of  surface  runoff  and 
baseflow  appear  by  time  periods  in  Table  4.  Figures  8 
and  9  illustrate  trends  for  the  estimated  percentages  of 
surface  runoff  and  baseflow  during  the  three  time 
intervals  of  the  study  period.  The  results  indicate  the 
relative  fraction  of  total  runoff  that  is  surface  runoff 
has  decreased  from  the  1940s  through  1980s.  Baseflow 
shows  an  opposite  trend.  Table  6  summarizes  results 
from  regression  analyses  performed  on  estimated  surface 
runoff  and  baseflow  percentages  during  the  study  period 
and  within  the  three  time  periods.  The  surface  runoff 
regression  line  for  the  study  period  is  Y  =  -0.0064X  + 
13.34;  baseflow:  Y=  0.0064X  -  12.34.  These 
statistical  relationships  are  significant  at  the  0.0000 
probability  level  (see  Table  6) .  The  variances 
explained  by  the  two  equations  are  low  (r2  =  0.24), 
which  is  evident  by  visual  inspection  of  Figures  8  and 
9 .  Separate  regression  am  yses  for  the  three  time 
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TABLE  6  -  SURFACE  RUNOFF ,  8ASEF19V,  AND  TOTAL  RUNOFF  VERSUS  TIRE 
FRON  1940  -  1987 


CORRELATION  OF  THE  PERCENT  OF  SURFACE  RUNOFF  VS  TIME:  -0.4811 
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CORRELATION  OF  THE  TOTAL  RUNOFF  VS  TINE:  -0.2841 
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TABLE  6  RE5JESSICN  ANALYSIS  -  SURFACE  RUNOFF  ANO  8ASEFL0D  VERSUS  TIME 
FOOD  i940  -  195! 
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TABLE  i  ’EGRESSION  ANALYSIS  -  SURFACE  RUNOFF  AKO  BASEFLBU  VERSUS  'INE 
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0.5012 

-0.3940 

-0  3435 

1.0000 

0.1716 

0.2U2 

-0.0915 

0.7325 

7  OAT  ART  PRECIP 

*0  0403 

-0.239! 

-0.2830 

-0.3210 

0.1716 

1.0000 

0.6360 

0  1 1*2 

0.1908 

U  OAY  ANT  PRECIP 

0.1S05 

-0 . 0543 

-0.2571 

-0.2S43 

0.2112 

0  6860 

1 .0000 

-0.0567 

0.2833 

RATIO  PRES/OUR 

-0.6930 

-0.0233 

-0.0373 

-0.2794 

-0.0915 

0.1112 

-OIS67 

1.0000 

0  0441 

J  SURFACE  RUNOFF 

0.0522 

0 . 4935 

-0  4279 

-0.4504 

0  7325 

0.1908 

0.333 

0.0441 

1.0000 

SE5HESSI0N  AN ALTS IS 

OErENOENT  VAR i ABLE’ 

PERCENT 

SURFACE  RUNOFF 

IROEPENOENT 

PASANE3R  STD: 3D 

STDflRD 

t-VALUE 

PR08ILTY 

SINPlE 

PARTIAL 

VARIABLES: 

ESTINATE  ESTINATE 

ERROR 

U*0) 

LEVEL 

R  SiR 

R  SOR 

INTERCEPT 

0:52'  0.0000 

0. !  *55 

2.6300 

0.0165 

STORN  DURATION 

-il  -y  4037 

0  1054 

-1 .6500 

0.1034 

0  .034 

0.  i29? 

STORN  :*£•:;? 

0  U'7.39  0  6014 

0  ';2?2 

2.5309 

0.0203 

0  2435 

0.2523 

1Y05Q  jU’fcTIOM 

tv.  00:0  0.140! 

0  -;0E4 

1  4700 

0 . 1 S66 

0  831 

0  1027 

TIME  TO  :>A* 

-ii  '■)  P2  -y  65:0 

0  007: 

-2.- SOO 

9. "447 

o ;;:? 

0  1355 

PEA-  OivCHAAGE 

0  VU2U  0  5'02 

0.00? 

2.U100 

0  0584 

0.5:56 

0  1753 

t  day  art  -re::; 

:>  '161  0  53? 

0  '.'iO'-J 

0,5000 

i)  4333 

o  '.Hi 

0  0326 

'4  ;ii  ART  -’EC.’ 

0.  .066  ■'  v?': 

i)  124 

.i  =400 

0  5S7' 

1  ..;m‘ 

0.914? 

RATIJ  i-ii  OUR 

-(■  S3?  ->  ?■?:•» 

•;  ;  ?01 

-i  4500 

0  '627 

,■  .3 

:» .  •? 
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«61E  6  REGRESSION  ANALYSIS  -  SURFACE  RUNOFF  ANO  BASEFLQU  VERSUS  TIRE 
FROil  1969  -  1987 


CORRaSTlDN  if  THE  PERCENT  OF  SURFACE  RUNOFF  VS  TIKE.  -0.0405 


STANOARD  t -VALVE  PROB  F 

-  RATIO 

PROB 

R 

REGRESSION  EQUATION: 

ERROR  U  =  0)  LEva 

LEVEL 

SQUARED 

PERCENT  OF  SURFACE  RUNOFF:  Y  =  -I.OOIS  -  S.S5 

0.0072  -0.2  0. Bill 

0.04 

0  Si! 

0.0016 

RECESSION  ANALYSIS  OF  INDEPENDENT  VARIA8LES 
FRCS  362  -  1968 


REGRESSION  ANALYSIS  F  -  RATIO  PROBABILITY  R  SQUARED 

SURFACE  RUNOFF  VS  INDEPENDENT  VARIABLES:  LEVEL 

2.8  4.035  9.5637 


CORRELATIONS 


STORN 

STORN 

HYDRO 

TIRE  TO 

PEAK 

7  OAY 

14  OAY 

RATIO 

SSURFACE 

DURATION 

PRECIP 

DURATION  PEAK 

Q 

ANT  PRE 

ANT  PRE 

PREC/DUR  RUNOFF 

STORN  DURATION 

1.0000 

0.5117 

0.0983 

0.4415 

0.3096 

-0.0686 

0.0497 

-0.7638 

0.1573 

STORN  PRECIP 

0.5117 

1 .0000 

0.4462 

0.3096 

0.5314 

0.0234 

0.0912 

-0.3344 

0.5300 

HYDRO  DURATION 

0.0983 

0.1462 

1.0000 

9.2545 

0.0294 

-0.3993 

-0.52S9 

-0.1621 

-0.0210 

TINE  TO  PEAR 

0.4415 

0.3095 

0.2545 

1.0000 

0.0275 

0.0014 

-0.0514 

-0.16)0 

0.0835 

PEAK  DISCHARGE 

0.3096 

0.5SW 

0.0234 

0.0275 

1.0000 

-0.1149 

0.3037 

-0.2026 

0.670! 

7  DAY  ANT  PRECIP 

-0.06S6 

0.0234 

-0.3993 

0.0014 

-0.1149 

1 .0000 

0.4153 

0.3269 

0.0996 

U  DAY  ANT  PRECIP 

0.0497 

0.0012 

-0.5289 

-0.0514 

0.3037 

0.4153 

1.0000 

0.0861 

0.2827 

RATIO  PREC/OUR 

-0.7603 

-0.3344 

-0.1624 

-0.1610 

-0.2026 

0.3269 

0.0861 

1.0000 

-0.1505 

t  SURFACE  RUNOFF 

0.1578 

0.5300 

-0.9210 

0.0835 

0.6701 

0.0996 

0.2827 

-0.1505 

1.0000 

REGRESSION  ANALYSIS 

OEfENOENT  VARIABLE: 

PERCENT  SURFACE  RUNOFF 

INDEPENDENT 

PARANETER 

STD I EEC 

STDARD 

t-VALVE 

PRQ81LTY 

SINPLE 

PARTIAL 

VARIABLES: 

ESTIMATE 

ESTIMATE  ERROR 

(1=9) 

LEVa 

R  SQR 

R  SQR 

INTERCEPT 

0.3403 

0.9909 

0. .754 

1.9400 

9.0691 

STOP"  DURATION 

-0.0074 

-0,4280 

0.0049 

-1.1900 

0.1512 

0.0219 

0.1157 

storn  =rk;p 

0  0607 

0.1127 

0.0429 

1.4200 

0.1744 

0.2309 

0.  DOST 

■>'!C:0  DURATION 

-v. 0019 

-0  709 

9.0023 

-0.6700 

0.5092 

0.0006 

0.0260 

TINE  '0  ?EAK 

0.0913 

0.9413 

0.006! 

0.2100 

9.3336 

0.3070 

9.9027 

'Em1  DISCHARGE 

o.oo s; 

0.5116 

0 . v0*3 

2. 1500 

0,0461 

0.1190 

0.2195 

7  DAi  -NT  :REC!P 

v.ois: 

0  -349 

0  0251 

0 . 6000 

0.5539 

0.D099 

0  0205 

.4  In.  ANT  VPECDP 

9  Ool? 

0 .9233 

0.0114 

0.1290 

0.3970 

0.0799 

0  0010 

RATIO  ’RED. -'UR 

-i  .is; 

-9  .90S 

0.1139 

-1.0230 

0.322S 

0.0226 

0  0576 
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periods  indicate  the  sharpest  decline  in  surface  runoff 
occurred  from  1940  to  1952  with  moderate  declines  in  the 
later  two  groups  (Table  6  and  Figure  8) .  The 
statistical  relationships  and  small  explained  variances 
for  all  groups,  especially  the  later  two,  are  not 
significant.  Nevertheless,  all  analyses  suggest  a 
significant  decrease  in  surface  runoff  from  the  1940s  to 
the  1980s. 

MULTIPLE  REGRESSION  ANALYSIS 

Multiple  regression  analysis  was  used  to  explore 
long-term  historical  trends  in  surface  runoff  and 
baseflow  with  trends  in  other  hydrologic  and  climatic 
factors  (Table  6) .  An  analysis  of  the  average  monthly 
precipitation  and  temperature  form  April  through  October 
for  each  year  of  the  study  period  indicates  wide  and 
random  fluctuations  (precipitation  r2=0.003, 
temperature  r2=0.019)  (Figure  10  and  11).  Analyses 
showed  no  statistical  significance  and  therefore  no 
trends  in  the  slopes  of  the  regression  lines  for  both 
precipitation  and  temperature.  The  monthly  average 
precipitation  data  from  April  through  October  were 
derived  from  the  previously  discussed  four  weather 
stations  using  the  Thiessen  method.  Temperature  data 
are  from  the  Darlington  station  which  has  the  only 


AVERAGE  WATERSHED  PRECIPITATION 


YEARLY  WATERSHED  TEMPERATURE 
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continuous  temperature  record  for  the  study  period.  The 
lack  of  change  in  precipitation  and  temperature  with 
time  indicate  that  climatic  factors  are  not  responsible 
for  the  observed  historical  trend  in  surface  runoff  and 
baseflow. 

Multiple  regression  analyses  were  used  for 
evaluating  climatic  influences  on  runoff.  Storm  and 
runoff  characteristics  constitute  the  independent 
variables  analyzed  in  relation  to  the  estimated  surface 
runoff  dependent  variable.  The  independent  variables 
are:  1.  storm  precipitation,  2.  storm  duration,  3. 

estimated  surface  runoff  duration,  4.  the  time  from  the 
beginning  of  surface  runoff  to  the  peak  discharge,  5. 
peak  discharge,  6.  preceding  seven  days  antecedent 
precipitation,  7.  preceding  fourteen  days  antecedent 
precipitation,  and  8.  storm  intensity  ratio  (used 
earlier  to  evaluate  similarity  between  storm  groups) 
(Table  4) . 

Durations  of  surface  runoff  were  estimated  from 
gage  strip  charts  or  from  tables  of  gage  heights  by 
identifying  the  time  between  the  initial  point  of  rise 
and  the  asymtotic  stabilization  point  on  the  recession 
limb.  Time  to  peak  data  was  determined  by  measuring  the 
time  from  the  beginning  of  surface  runoff  to  the  time 
the  highest  peak  observed  on  the  hydrograph.  This 
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method  was  used  instead  of  measuring  from  the  centroid 
of  rainfall  to  the  centroid  of  the  hydrograph  due  to  the 
inaccuracy  of  the  rainfall  centroid.  The  rainfall 
centroid  is  inaccurate  due  to  the  use  of  the  Thiessen 
method  and  having  only  one  hourly  recording 
precipitation  station.  Many  hydrographs  have  multiple 
peaks.  To  ensure  consistency  the  highest  magnitude  peak 
was  used  to  determine  estimates  of  the  time  to  peak  and 
peak  discharge  variables.  Use  of  the  highest  magnitude 
hydrograph  peak  is  desired  for  consistency  and 
introduced  slight  fluctuations  in  the  time  to  peak 
data. 

The  Thiessen  method  was  used  to  determine  both 
antecedent  precipitation  variables.  The  definition  of 
the  seven-day  antecedent  precipitation  was  based  on 
recommendations  given  in  the  National  Engineering 
Handbook  and  by  hydrologists  of  the  U.S.  Soil 
Conservation  Service  (S.C.S.)  (Mockus,  1972  and  phone 
conservation  with  John  Milligan,  S.C.S. ,  Madison,  Wl, 
June  1989) .  The  National  Engineering  Handbook 
establishes  a  criteria  of  five  days  of  no  antecedent 
precipitation  to  achieve  maximum  infiltration  capacity 
while  hydrologists  at  the  Wisconsin  S.C.S.  believe  at 
least  seven  days  are  needed.  Following  the  same 
procedures  and  logic,  a  fourteen-day  antecedent 
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precipitation  criterion  was  used  to  evaluate  longer-term 
influences . 

Table  6  summarizes  the  results  of  the  multiple 
regression  analysis.  Results  indicate  weak  correlations 
between  other  climatic  and  hydrologic  variables  and  the 
percentage  of  estimated  surface  runoff  except  for  peak 
discharge.  The  strong  relationship  between  peak 
discharge  and  the  percent  of  estimated  surface  runoff 
(r=0.70)  indicates  a  simultaneous  decline  in  peak 
discharge  in  association  with  declining  magnitude  for 
the  fraction  of  total  runoff  that  is  surface  runoff. 

Time  trends  for  average  peak  discharges  support  the 
results  from  multiple  regression  analyses  (Table  4) . 

The  strongest  correlation  coefficient  (r=0.75)  occurred 
between  storm  duration  and  the  dimensionless  ratio  of 
storm  precipitation  verses  storm  duration.  The 
relationship  is  not  causally  significant  because  storm 
duration  is  a  major  component  of  the  ratio.  Similarly, 
the  strong  relationship  between  the  seven-day  and 
fourteen-day  antecedent  precipitation  (r=0.61)  occurs 
because  the  seven-day  precipitation  is  part  of  the 
fourteen-day  precipitation. 

Regression  analysis  showed  the  equation  linking  the 
independent  variables  and  the  percent  of  estimated 
surface  runoff  is  statistically  significant  to  the 
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0.0000  level  with  59%  of  the  variance  explained  (Table 
6) .  The  majority  of  the  variance  is  attributed  to  peak 
discharge  (partial  r2  =  0.36).  The  parameter 
estimates  have  a  relative  high  standard  error  for  all 
independent  variables.  Statistically  significant 
variables  are  peak  discharge  (0.0000),  time  to  peak 
(0.0037),  and  precipitation  amount  (0.0019).  The  strong 
positive  association  between  peak  discharge  and  percent 
estimated  surface  runoff  indicates  that  both  have  been 
decreasing  since  the  1940s.  The  statistical 
significance  between  the  time  to  peak  and  percent 
estimated  surface  runoff  may  be  questionable  due  to  the 
influence  of  multiple  peaked  hydrographs  as  discussed 
earlier.  Nevertheless,  the  weak  inverse  relationship  is 
consistent  with  the  notion  that  surface  runoff  would 
likely  be  greatest  for  storms  with  flashy  hydrographs 
with  short  lag  times  between  rainfall  and  runoff.  The 
relationship  between  precipitation  amount  and  percent 
estimated  surface  runoff  is  weaker  than  what  might  be 
expected  although  the  relationship  does  demonstrate  that 
the  proportion  of  total  runoff  that  is  surface  runoff 
tends  to  increase  as  storm  magnitude  increases. 
Statistical  analyses  indicates  that  climatic  differences 
occurring  between  the  1940s  and  the  1980s  are  not 
responsible  for  the  decreasing  trend  in  surface 
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runoff. 

A  regression  analysis  was  used  to  link  specific 
climatic  and  land  use  variables  to  significant 
hydrologic  results.  Limiting  the  number  of  independent 
variables  decreased  the  influence  from  a  large  number  of 
independent  variables  and  allowed  a  clearer 
interpretation  of  specific  climatic  and  land  use 
influences.  Climatic  variables  included  storm 
precipitation,  storm  duration,  and  seven-day  antecedent 
precipitation.  The  land  use  variable  used  was  the 
yearly  percent  of  land  under  soil  conservation 
practices.  Yearly  percentages  of  land  under  soil 
conservation  practices  were  derived  from  the  seven 
data  sets  observed  from  aerial  photographs  (see  Table 
5) .  The  yearly  percentages  was  assumed  to  vary  linearly 
between  aerial  photograph  data  sets.  Significant 
hydrologic  results  included  estimated  percent  of  surface 
runoff,  time  to  peak,  and  peak  discharge. 

Regression  analyses  indicate  strong  relationships 
between  the  independent  variables  of  storm  precipitation 
and  land  use  and  the  dependent  variables  of  peak 
discharge  and  estimated  percent  of  surface  runoff  (see 
Table  7) .  The  equation  for  the  peak  discharge  dependent 
variable  was  significant  at  the  0.000  level  with  storm 
precipitation  and  land  use  significant  at  the  0.000 
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TABLE 


RESRES 


sfATERSHEO  AND  STORM  ASPECTS  1940  -  1987 


TO  PEA!::  VS  INDEPENDENT  VARIABLES 


REGRESSION  ANALYSIS 

F  -  RATIC 

1  PROBABILITY  R 

SQUARED 

TINE  TO  PEAK  VS  INDEPENDENT  VARIABLES: 

LEVEL 

2.82 

0.041 

0.1185 

CORRELATIONS 

STORM 

STORM  7  DAY 

X  LAND  US1N6 

TIME  TO 

DURATION 

PRECIP  ANT  PRE 

SOIL  CONSERVATION 

PEAK 

STORN  DURATION  1.0000 

0.4075  -0.0928 

-0.0496 

0.2866 

STORM  PRECtP  0.4075 

1.0000  -0. 1288 

0.0594 

0.2587 

7  DAY  ANT  PRECIP  -0.0928 

-0.1268  1.0000 

-0.0951 

-0.1345 

I  LAND  US I NS 

SOIL  CONSERVATION  -0.01% 

0.0596  -0.0951 

1.0000 

0.0661 

TIME  TO  PEAK  0.2866 

0.2587  -0.1345 

0.0661 

1.0000 

REGRESSION  ANALYSIS 

DEPENDENT  VARIABLE:  TIME  TO  PEAK 

INDEPENDENT 

PARAMETER  STOIZED 

STDARD  t-VALUE 

PROBILTY 

SIMPLE 

PARTIAL 

VARIABLES: 

ESTIMATE 

ESTIMATE  ERROR  (1=0) 

LEVEL 

R  SfJR 

R  SQR 

INTERCEPT 

5.4558 

0.0000 

3.0448  1.7900 

0.0770 

STORM  DURATION 

0.1574 

0.2182 

0.0844  1.8700 

0.0659 

0.0822 

0.0427 

STORM  PRECIP 

0.3450 

0.1550 

0.6396  1.3200 

0.1903 

0.0669 

0.0219 

7  DAY  ANT  PRECIP 

-0.3363 

-0.0890 

0.4071  -0.8300 

0.4114 

0.0181 

0.0087 

X  LAND  USING 

SOIL  CONSERVATION 

0.0343 

0.0S93 

0.0621  0.5500 

0.5824 

0.0044 

0.0039 

TABLE  7  REGRESSION  ANALYSIS  -  WATERSHED  AND  STORM  ASPECTS  1940  -  1987 
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AK  DISCHARGE  VS  INDEPENDENT  VARIABLES 


REGRESSION  ANALYSIS  F  -  RATIO 

PEAK  DISCHARGE  VS  INDEPENDENT  VARIABLES: 

9.05 

CORRELATIONS 


STORM 

STORM 

7  DAY 

%  LAND  USING 

PEAK 

DURATION 

PRECIP 

.  ANT  PRE 

SOIL  CONSERVATION 

OISCHARGE 

STORM  DURATION  ! .0000 

0.4075 

-0.0928 

-0.0496 

0.1495 

STORM  PRECIP  0.407S 

1.0000 

-0.1268 

0.0594 

0.3584 

7  DAY  ANT  PRECIP  -0.0923 

-0.1268 

! .0000 

-0.0951 

0.2245 

S  LAND  USING 

SOIL  CONSERVATION  -0.0496 

0.0596 

-0.0951 

1.0000 

-0.3390 

PEAK  DISCHARGE  0.1 49S 

0.3584 

0.2245 

-0.3390 

1.0000 

PROBABILITY  R  SQUARED 
LEVEL 

0.000  0.317! 


REGRESSION  ANALYSIS 
DEPENDENT  VARIABLE:  PEAK  DISCHARGE 


INDEPENDENT 

PARAMETER 

STDIZEO 

STOARD 

t-VALUE 

PROBILTY 

SIMPLE 

PARTIAL 

VARIABLES: 

ESTIMATE 

ESTIMATE  ERROR 

<X=0> 

LEVEL 

R  SQR 

ft  SQR 

INTERCEPT 

-14.6010 

0.0000 

13.1761 

-1.1100 

0.2712 

STORM  DURATION 

-0.0501 

-0.0141 

0.3651 

-0.1400 

0.8913 

0.0223 

0.0002 

STORM  PRECIP 

11.1286 

0.4)53 

2.7676 

4.0200 

0.0001 

0.1285 

0.1717 

7  DAY  ANT  PRECIP 
l  LAND  USING 

4.5228 

0.2434 

1.7620 

2.5700 

0.0122 

0.0504 

0.0779 

SOIL  CONSERVATION 

-0.9715 

-0.3412 

0.2689 

-3.6100 

0.0005 

0.1149 

0.1433 
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TABLE  7  REGRESSION  ANAL VS IS  -  WATERSHED  AND  STORM  ASPECTS  1340  -  1387 


ORFACE  RiJNSFF  VS  INDEPENDENT  VARIABLES 

REGRESSION  ANALYSIS  F  -  RATIO 

X  SURFACE  RUNOFF  VS  INDEPENDENT  VARIABLES: 

17.40 


uORRtLAi iGNS 


PROBABILITY  R  SQUARED 
LEVEL 

0.000  0.4715 


STORM 

STORM 

7  DAY 

X  LANO  US1N6 

*  SURFACE 

DURATION 

PRECIP 

ANT  PRE  SOIL  CONSERVATION 

RUNOFF 

STORM  DURATION 

1  0000 

0.4075 

-0.0328 

-0.0436 

0.1435 

STORM  PRECIP 

0 . 4075 

1.0000 

-0.1268 

0.0534 

0.3584 

7  DAY  ANT  PRECIP 

-0  0328 

-0.1268 

1.0000 

-0.0351 

0.2245 

5  LAND  OSINS 

SOIL  CONSERVATION 

-0.04S8 

0.0536 

-0.0351 

1.0000 

-0.3330 

X  SURFACE  RUNOFF 

0.0474 

0.3542 

0.2113 

-0.5087 

1.0000 

REGRESSION  ANALYSIS 
DEPENDENT  VARIABLE:  1!  SURFACE  RUNOFF 


INDEPENDENT 

PARAMETER 

STDIZED 

STDARO 

i-VALUE 

PROBILTY 

SIMPLE 

PARTIAL 

VARIABLES:' 

.ESTIMATE 

ESTIMATE  ERROR 

(7=0) 

LEVEL 

R  SQR 

R  SQR 

INTERCEPT 

0.4256 

0.0000 

0.0683 

6.2300 

0.0000 

STORM  DURATION 

-0.0032 

-0.1524 

0.0013 

-1.6800 

0.0964 

0.0023 

0.0350 

STORM  PRECIP 

0.0784 

0.4738 

0.0144 

5.2100 

0.0000 

0.12S4 

0.2585 

7  DAY  ANT  PRECIP 

S  LANO  USING 

0.0228 

0.2079 

0.0091 

2.4900 

0.0148 

0.0443 

0.0738 

SOIL  CONSERVATION 

-0.0088 

-0.5246 

0.0014 

-6.3100 

0.0000 

0.2588 

0.3383 

74 


level.  Storm  precipitation  and  land  use  had  the  highest 
correlations  with  peak  discharge.  The  equation  for  the 
estimated  percent  of  surface  runoff  was  significant  at 
the  0.0000  level  and  had  relationships  with  storm 
precipitation  and  land  use  similar  to  peak  discharge. 
Indications  from  regression  analyses  support  the  link 
between  land  use  and  the  general  decreasing  trends 
observed  in  surface  runoff  and  peak  discharge.  Storm 
precipitation  appears  to  fluctuate  with  the  percent  of 
surface  runoff  and  peak  discharge  as  expected. 

Regression  results  indicated  no  statistical  significance 
between  the  independent  variables  and  time  to  peak 
dependent  variable.  The  equations  explained  above  could 
be  reformulated  eliminating  the  insignificant 
independent  variables  which  would  allow  a  better 
interpretation  of  the  significant  variables.  The  above 
process  was  done  to  illustrate  the  insignificant 
influence  from  variables  likely  to  influence  the 
dependent  variables. 

STORM  GROUP  COMPARISON 

The  similarity  of  climatic  conditions  throughout 
the  study  period  from  the  1940s  through  the  late  1980s 
provides  a  framework  for  evaluating  the  influence  that 
soil  conservation  practices  have  had  on  regulating  the 
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relative  fractions  of  surface  runoff  and  baseflow 
runoff. 

The  percentage  of  total  runoff  that  was  delivered 
as  surface  runoff  experienced  more  than  a  20%  decrease 
between  the  1940s  and  1980s  (Table  4  and  Figure  8) . 

Much  of  the  reduction  may  have  occurred  between  the 
1960s  and  1980s  period.  Separate  regression  equations 
describing  trends  of  surface  runoff  magnitude  within  the 
three  groups  are  not  statistical  significant  due  to  wide 
variance' within  the  groups.  However,  Figure  8 
illustrates  the  decreasing  trend  in  surface  runoff 
within  the  three  periods. 

Figure  8  data  suggest  a  sharp  reduction  in  the 
percentage  of  estimated  surface  runoff  between  the  1960s 
and  1980s.  This  reduction  coincides  with  an  equally 
sharp  increase  in  the  observed  percentages  of  land  area 
using  soil  conservation  practices  (Table  5  and  Figure 
7) .  However,  it  is  clear  that  the  gradual  reduction  in 
the  fraction  of  total  runoff  that  is  surface  runoff 
between  the  1940s  and  1980s  coincides  with  a  gradual 
increase  in  soil  conservation  practices. 

An  evaluation  of  responses  of  surface  runoff  to 
specific  ranges  of  rainfall  magnitude  showed  a 
decreasing  trend  within  all  ranges  except  within  the 
greatest  rainfalls  (Figures  12  -  14) .  The  continued 


ire  12.  The  estimated  percent  of  total  runoff  that 
surface  runoff  from  storms  that  had  a  precipitatlo 
int  between  2.t>0  and  4.00  centimeters. 
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high  fraction  of  total  runoff  that  is  surface  runoff 
occurring  within  the  6.00-8.00  cm  rainfall  (Figure  14) 
is  expected.  High  percentages  of  surface  runoff  occur 
when  rainfall  contributes  exclusively  to  surface  runoff 
after  the  storage  capacity  of  the  land  is  exceeded.  The 
outliers  in  Figure  14  probably  occurred  as  a  result  of 
dry  antecedent  conditions  (1984  outlier)  or  the 
influence  of  localized  intense  convectional  storms  to  be 
under  represented  by  the  Thiessen  method  (1947 
outlier) .  There  were  no  apparent  causes  for  the  1974 
outlier. 

The  time  to  peak  independent  variable  increases 
among  the  three  time  periods  supporting  the  general 
trend  in  surface  runoff  (Figure  15) ,  but  analyses  of 
variances  among  the  three  groups  indicate  no  statistical 
significance  (Table  8) .  An  insignificant  relationship 
is  expected  because  watershed  physiography  exerts  the 
dominant  role  on  surface  runoff  during  this  study. 

The  average  peak  discharge  decreased  64%  during  the 
study  period  (Table  4) .  Analysis  of  variance  indicates 
statistical  significance  to  the  0.01  level  among  the 
three  time  periods,  to  the  0.03  level  between  the  1960s 
and  1980s,  and  no  statistical  significance  (0.19 
probability  level)  between  the  1940s  and  1960s 
supporting  the  previous  indication  of  similarity  between 
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TABLE  5 

»‘*C 

iv.J 

OF  VARIANCE 

STANDARD  ■ 

PERIOD 

DEVIATION  F 

-RATIO 

PR03A8ILITY 

MEANS 

TIN£  *3  PEAK ' 

ALi.  THREE  PERIODS 

* :  jc 

•j-j 

6.30 

»  5 

V  .  **J 

0.62 

PEilOil  !  V‘:>  i. 

0.44 

0.51 

PERIOD  5; 

l  7  *. 

PEE i CD  1  *3  3ER!OD  3 

0  1 

0 . 75 

PERIOD  2: 

I'i  .34 

:cC'T,',r,  i  t;c  occtop  ? 
itWdu  t  f  j  •  ur.t'ju  y 

0  39 

0.35 

PERIOD  3: 

12.33 

-EAK  DISCHARGE: 

ALL  THREE  PERIODS 

iv  .  v8 

29  S-F 

5.15 

0 . 008 

PERIOD  !  VS  PERIOD  0 

1.69 

0.19 

PERIOD  l; 

40.67 

PERIOD  C  VS  PERIOD  3 

4  ?? 

0.03 

PERIOD  2: 

28.77 

PERIOD  i  VS  PERIOD  S 

1 A  I  C 

1  W  .  13 

0.002 

PERIOD  3! 

!  4  80 

7  CAY  ANTECEDENT 

PRECIPITATION: 

i  .52 

1  .69 

0  47 

0.63 

4  DAY  ANTECEDENT 

PRECIPITATION: 

3.37 

2.47 

0.29 

0.75 

HY0RQ6RAPH  DURATIONS: 

ALL  THREE  PERIODS 

45.26 

IS. 12 

1.35 

0.164 

PERIOD  1  VS  PERIOD  2 

!  .04 

0.31 

PERIOD  l : 

45.06 

PERIOD  2  VS  PERIOD  3 

3.99 

0.051 

PERIOD  2: 

49.24 

PERIOD  !  VS  PERIOD  3 

0.72 

0.4 

PERIOD  3: 

41.48 

PERCENTAGE  OF 

SURFACE  RUNOFF: 

all  three  periods 

0.62 

0.16 

16.07 

0.000 

PERIOD  i  v'S  PERIOD  2 

0.5 

0.484 

PERIOD  l: 

71* 

PERIOD  2  VS  PERIOD  3 

19.31 

0.000 

PERIOD  2: 

68* 

PERIOD  1  VS  PERIOD  3 

27.35 

0.000 

PER I 00  3: 

48* 

PERCENTAGE  OF  8ASEFL0V: 

ALL  three  periods 

0  36 

0.16 

16.07 

o.ooo 

PERIOD  1  VS  PER I DO  2 

0.5 

0.484 

PERIOD  l: 

29* 

PERIOD  2  VS  PERIOD  3  : 

- 

19.3! 

0.000 

PERIOD  2: 

32  S 

PERIOD  1  VS  PERIOD  3 

27.35 

0.000 

PERIOD  3: 

52S 
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the  two  earlier  time  periods  (Table  8) .  The  average 
seven-day  and  fourteen-day  antecedent  precipitation 
among  the  three  time  periods  are  not  significantly 
different  statistically  (Table  8) . 

The  average  durations  of  surface  runoff  among  the 
three  time  periods  were  similar  with  no  statistical 
significance  indicated  at  the  0.05  level  (see  Tables  4 
and  8).  However,  a  statistically  significant  (0.10 
level)  difference  occurs  between  the  1960s  and  1980s. 
The  overall  similarity  in  durations  of  surface  runoff 
among  the  three  time  periods  coupled  with  similarities 
in  storm  intensities  should  produce  similar  hydrographs 
according  to  the  unit  hydrograph  theory  if  watershed 
environmental  conditions  are  constant  (Dunne  and 
Leopold,  1978) .  Figure  15  shows  composite  unit 
hydrographs  for  each  time  period  generated  from  average 
group  hydrograph  characteristics.  Note  that  peak 
discharge  decreases  by  64%  from  the  1940s  to  the  1980s, 
while  the  time  to  peak  discharge  correspondingly 
increases  from  10.72  hours  in  the  1940s  to  12.38  hours 
in  the  1980s. 

DISCUSSION 

The  changing  characteristics  of  the  composite  unit 
hydrographs  in  Figure  15  reflect  the  increased 
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application  of  soil  conservation  practices  between  the 
1940s  and  1980s.  The  increased  use  of  soil  conservation 
practices  are  responsible  for  the  smaller  peak  flows  and 
longer  durations  for  time  to  peak  discharges.  The 
decrease  in  peak  discharge  and  the  increase  in  time  to 
peak  discharge  further  support  the  view  that  less  of  the 
total  runoff  from  precipitation  in  the  1980s  is  being 
allocated  to  surface  runoff  than  was  characteristic  of 
the  1940s.  Although  by  the  method  of  hydrograph 
separation  employed  here,  baseflow  is  a  determined 
complement  of  surface  runoff,  it  is  useful  to  note  that 
average  baseflow  went  from  25%  of  total  runoff  in  the 
1940s  to  45%  of  total  runoff  in  the  1980s. 

The  results  of  the  present  research  show  that 
watershed  precipitation  is  being  retained  and  released 
more  slowly  in  the  1980s  than  it  was  in  the  1940s,  and 
that  the  difference  is  due  to  the  increased  utilization 
of  soil  conservation  practices.  A  plot  of  total  runoff 
from  storm  events  indicates  a  gradual  decrease  in  total 
runoff  during  the  study  period  (Figure  16  and  Table  6) . 
The  decrease  in  total  runoff  occurs  while  storm 
precipitation  remains  constant  (Figure  6)  further 
indicating  that  improved  soil  conservation  land  use  has 
increased  retention  of  precipitation  within  the 
watershed  during  the  later  half  of  the  study  period.  An 


Figure  16.  The  total  amount  of  runoff  which  occurred 
from  the  83  analyzed  storms.  The  output  Is  cubic  meters 
as  discussed  In  Chapter  Five. 


apparent  similarity  exists  between  the  first  and  second 
time  periods  average  percent  of  estimated  surface 
runoff. 

CONCLUSIONS 

The  results  of  this  study  reinforce  conclusions  of 
previous  studies  that  examined  hydrologic  response  co 
land  use  changes  within  the  Driftless  area  (Knox,  1972; 
1977;  1987;  Trimble  and  Lund,  1982;  and  Sartz  1970). 
These  studies  showed  that  major  reductions  in 
sedimentation  and  surface  runoff  occurred  when  soil 
conservation  practices  were  introduced  primarily  after 
1950. 

Most  previous  studies  have  estimated  hydrologic 
responses  to  land  use  and  conservation  practices  by 
observing  historical  changes  in  sedimentation,  channel 
morphology,  or  by  measuring  runoff  from  smalt  plots 
under  varying  land  use.  This  study  has  focused  on 
quantitative  evaluations  of  streamflow  hydrographs  to 
measure  the  fractions  of  surface  and  baseflow  runoff 
from  a  medium  sized  basin. 

The  results  of  the  present  study  show  that  the 
fraction  of  total  runoff  that  is  delivered  as  surface 
runoff  decreased  from  the  1940s  to  the  1980s  in  the 
Galena  watershed.  Statistical  analyses  of  climatic  data 
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showed  that  the  reduction  was  not  related  to  climatic 
causes.  Furthermore /  physiographic  influences  other 
than  land  use  also  remained  constant  over  the  short 
period  of  this  study.  On  the  other  hand,  land  use 
changed  significantly,  especially  in  the  adoption  of 
soil  conservation  practices  within  the  Galena  watershed 
between  the  1940s  and  the  1980s.  Quantitative 
comparison  of  temporal  changes  in  soil  conservation 
practices  and  runoff  hydrology  showed  that  increased 
percentages  of  land  under  soil  conservation  practices  is 
mainly  responsible  for  the  observed  decrease  in  surface 
runoff  and  flood  peaks. 
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CHAPTER  6 :  SUMMARY 

Land  use  within  the  Galena  watershed  from  1940 
through  1987  changed  as  farmers  expanded  their  use  of 
soil  conservation  practices.  The  repercussion  from 
increased  soil  conservation  has  influenced  the  various 
fractions  of  total  runoff  decreasing  surface  runoff 
while  increasing  baseflow  and  retaining  moisture  longer 
within  the  watershed. 

Prior  to  the  1940s  land  use  throughout  southwestern 
Wisconsin  suffered  from  extensive  soil  erosion  resulting 
from  traditional  agricultural  practices.  Beginning  in 
the  mid-1930s  soil  conservation  began  to  appear  and 
traditional  agricultural  practices  changed  as  a  result 
of  educational  efforts  by  the  Soil  Conservation  Service 
and  the  University  of  Wisconsin.  As  soil  conservation 
increased,  the  rate  of  soil  erosion  subsided  and  the 
surface  hydrology  changed.  Many  previous  studies 
discuss  influences  that  agricultural  practices  had  on 
surface  hydrology  within  the  Driftless  area  by  comparing 
sedimentation  rates,  channel  morphology,  and  estimating 
runoff  between  the  pre  and  post  1950s.  This  study 
attempts  to  focus  on  hydrologic  changes  occurring  from 
1940  through  1987  by  estimating  from  hydrographs  the 
percentages  of  total  runoff  that  is  surface  runoff  and 
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basef low. 

Storms  within  a  specified  range  were  selected  for 
analyses  to  accurately  evaluate  the  influence  from  soil 
conservation.  Land  use  analyses  defined  three  periods 
of  varying  soil  conservation  activity.  Statistical 
analyses  indicated  that  storm  intensities,  durations, 
and  magnitudes  were  not  statistically  different  among 
the  three  time  periods  and  throughout  the  study  period 
for  the  controlled  data  of  this  study.  A  slight  change 
in  storm  duration  was  observed  when  all  data  were 
included.  The  use  of  the  controlled  data  eliminated 
climate  as  a  cause  of  any  hydrologic  change  in  the 
surface  runoff  proportion  of  total  runoff. 

Land  use  was  evaluated  by  analyzing  seven  sets  of 
aerial  photographs  taken  in  1937,  1951,  1955,  1962, 

1968,  1976,  and  1985.  Aerial  photographs  are  of 
sufficient  quality  to  allow  land  use  classifications  of 
cultivated,  non-cultivated,  forested,  and  the  percent  of 
cultivated  land  under  soil  conservation.  The  scales  of 
the  aerial  photographs  did  not  allow  classification  of 
wide  row  verses  close  row  crops. 

Soil  conservation  encompasses  many  aspects.  The 
aspects  focused  on  in  this  study  are  primarily  contour 
farming  and  strip  cropping.  Terracing  and  the  use  of 
buffer  strips  were  difficult  to  differentiate  from 


contour  fanning  and  strip  cropping  and  were  included  as 
soil  conservation  adoption.  The  percent  of  cultivated 
and  non-cultivated  land  remained  nearly  constant  over 
the  study  period.  Small  fluctuations  occurred  in  the 
percent  of  forested  land.  The  most  significant  change 
occurred  in  the  increase  in  the  amount  of  land  using 
soil  conservation  practices. 

Percentages  of  total  runoff  that  was  surface  runoff 
and  baseflow  resulting  from  given  storms  were  determined 
from  hydrograph  analyses.  Surface  runoff  was  visually 
determined  from  storm  hydrographs  beginning  with  the 
initial  inflection  on  the  rising  hydrograph  limb  and 
terminating  with  the  asymtotic  stabilizing  point  on  the 
recession  limb.  A  data  file  defining  each  hydrograph  in 
(X, Y)  coordinates  was  then  processed  through  a  Fortran 
computer  program  which  determined  an  estimate  for  the 
percentages  of  total  runoff  that  was  surface  runoff  and 
baseflow.  The  resulting  percentages  served  as  dependent 
variables  in  equations  related  to  other  hydrologic  and 
climatic  independent  variables.  Multiple  regression 
analysis  and  analysis  of  variance  were  used  to  evaluate 
independent  variables  determined  from  storm  and 
hydrograph  characteristics. 

In  all  analyses  climatic  influences  appear  to  have 
weak  to  non-significant  relationships  with  the  fraction 
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of  total  runoff  that  was  surface  runoff  as  expected  due 
to  the  control  placed  on  storm  selection  and  relatively 
long-term  climatic  stability  within  the  Galena  watershed 
during  the  study  period.  Physiographic  factors 
influencing  runoff  except  for  land  use  remained  constant 
during  the  short  time  period  of  this  study.  Results  of 
the  analyses  showed  that  the  fraction  of  total  runoff 
that  is  surface  runoff  decreased  by  an  average  of  20% 
between  the  1940s  and  1980s.  A  necessary  corresponding 
increase  in  baseflow  therefore  also  resulted.  The  only 
significant  watershed  change  during  the  study  period 
that  could  explain  these  hydrologic  changes  was 
increased  use  of  soil  conservation  practices. 

An  accelerated  rate  of  adoption  of  soil 
conservation  practices  began  in  the  late  1960s  and  has 
continued  to  the  present.  A.S.C.S.  officials  estimate 
the  percent  of  cultivated  land  utilizing  soil 
conservation  may  be  as  high  as  95%  before  1995  due  to 
federal  legislation  defining  soil  conservation  programs 
as  prerequisites  to  receiving  federal  subsidies. 
Additional  research  conducted  in  the  mid-1990s  on  runoff 
within  the  Galena  watershed  should  be  undertaken  to 
document  the  hydrologic  response  to  this  anticipated 
land  use  change. 
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APPHTOIX  1 


FQRTRAJI  PROGRAM 


'ISO'S  J75 
TOHN  P  BAKER 
i A  MCV  SB 
LAB  #  A 


C  THE  STATELY  OF  THIS  PROGRAM  IS  TO  CREATE  A  COMMON  BLOCK  TO  STORE  DATA 
C  THAT  IS  ACCESSIBLE  TO  DETERMINE  THE  AREA  OF  IRREGULARILY  SHARED  POLr'- 
C  GCNS.  THIS  PROGRAM  WILL  ALSO  DETERMINE  BASIC  STATISTICS  IN  A  SUB- 
C  ROUTINE.  THIS  PROGRAM  WORKS  ON  AN  X  AND  Y  COORDINATE  SYSTEM.  TO 
C  FIND  THE  AREA  OF  YOUR  POLYGON,  THE  CORNERS  OF  YOUR  POLYGON  MUST  BE 
C  ON  AN  X  AND  Y  COORDINATE  SYSTEM.  STATISTICS  ARE  DETERMINED  FROM  A 
C  DATA  SET  YOU  INPUT  FOLLOW  THE  STATS  SUBROUTINE  FOR  THE  DATA  BASE 
C  FORMAT 

PROGRAM  POLYGON 

CHARACTER  S  I  0  NAME 

INTEGER  CHOICE.  INPUT.  COUNT 

REAL  0 1 GRAY  ( ! OOO , 5  > 

COMMON  D I GRAY,  COUNT 

COUNT**) 

PRINT*. 'TO  DETERMINE  THE  AREA  OF  AN  IRREGULARILY  SHAPED 
PRINT*, 'POLYGON,  TYPE  1.  TO  DETERMINE  STATISTICS  TYPE  0. 

PRINT*, 'TYPE  IN  YOUR  CHOICE  NOW: 

PRINT*,' 

READ'S, 100)  CHOICE 
100  FORMAT'!  14) 

IF  (CHOICE  E'5.  2)  THEN 
GO  TO  10 
END  IF 

600  PRINT*, ’BEFORE  BEGINNING,  ENSURE  YOUR  DATA  FILE  IS  ORGANISED  IN' 
PRINT*, 'THE  FOLLOWING  MANNER:  ' 

PRINT*, '  ' 

PRINT*,'  1.  A  DATA  FIELD  IDENTIFIER  PLACED  IN  THE  FIRST  ' 

PRINT*, 'FIVE  COLUMNS  ' 

PRINT*,’  ' 

PRINT*,'  2.  PLACE  A  FLAG  OF  ?•??■?  AT  THE  BEGINNING  AND  END' 

PRINT*,' OF  EACH  DATA  FIELD  IN  THE  FIRST  FIVE  COLUMNS  ' 

PRINT*, '  ' 

PRINT*,'  3.  YOUR  X  COORDINATE  IN  THE  NEXT  ELEVEN  COLUMNS  ' 

PRINT*, '  ' 

PRINT*,'  4.  YOUR  Y  COORDINATE  IN  THE  NEXT  ELEVEN  COLUMNS' 

PRINT*, 'FOLLOWING  THE  X  COORDINATE.' 

PRINT*, ’  ' 

PRINT*, 'IF  EVERYTHING  IS  READY,  PLEASE  ENTER  THE  NAME  OF  YOUR" 
PRINT* , ' DATA  FILE  (LIMIT  IT  TO  TWELVE  CHARACTERS). 

READ'S, 200)  NAME 
200  FORMAT' AID) 

WR I TE i 2 , *  )  ( NAME ) 

C  DETERMINING  SIDE  OF  INPUTTED  DATA  FILE 


OPEN' UNIT*),  FILE=NAME,  STATUS* ' OLD ' ) 
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1 =  READ I  1  .  300 .  END* 30) 

y.:o  format <  fs  o,  :fu 

CG'.'NT =CQUNT  1 
30  TO  1 5 

30  CONTINUE 
REWIND'  1  ) 

C  FILLING  ARRAY 

00  3-1  1  =  1.  COUNT 

READ'  1,300)  DIGRAYc  I .  1  > ,  OIGRAYt  I  ,  G  >  ,  01  GRAY <  1,3) 

31  CONTINUE 
RE'JINOi  1  ) 

CALL  AREA 
30  TO 

10  CALL  STATS 
GO  TO  333 

333  PRINT*, ’TYPE  1  IF  YOU  WANT  TO  DETERMINE  ANOTHER  AREA,  0  IF  YOU’ 
PRINT*. 'WANT  TO  DETERMINE  STATISTICS,  OR  3  TO  END  THIS  PROGRAM  ’ 
READ'S, 500)  INPUT 
500  FORMAT  CIA) 

IF  i INPUT  EQ.  1  )  THEN 
C0UNT=0 
GO  TO  600 

ELSE IF  C INPUT  . EQ .  2)  THEN 
GO  TO  10 
ELSE 

CONTINUE 
END  IF 
END 

SUBROUTINE  AREA 

REAL  DIGRAY  (1000,3),  AREAS.  AREA4 .  AREAS,  AREAS ,  AREAT ,  ARE ATI 
INTEGER  COUNT,  NCT ,  ID 
COMMON  01 GRAY,  COUNT 

PRINT*, 'WE  CAN  NOW  DETERMINE  THE  AREA  OF  YOUR  IPREGULARILY 
PRINT*, 'SHAPED  POLYGON.  ENSURE  YOUR  DATA  IS  NOT  LINKED  TO-’ 
PRINT*, 'GETHER  (THE  LAST  TERM  IS  NOT  THE  SAME  AS  THE  FIRST)  ' 
PRINT*, 'PLEASE  TYPE  IN  THE  OATA  IDENTIFICATION  NUMBER  ( THE  ' 
PRINT*, 'FIRST  FIVE  COLUMNS  OF  YOUR  DATA  FIELD)  NOW. 

READ'S, 101)  ID 
101  FORMAT (14) 

NCT*0 
AREA4=0  0 
AREA5-0  0 
AREA6=0  0 

C  NEED  TO  DETERMINE  THE  LAST  DATA  POINT  IN  THIS  FIELD  FOR  AREA 

00  10  1=1, COUNT 

IF  ( D I  GRAY  (I,l‘>  EO  ID)  THEN 
NCT=NCT  +  1 
END  IF 
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1  i.i  CONTINUE 

C  DETERMINING  AREA.  THE  FIRST  AND  LAST  ALGEBRAIC  TERMS  ARE  DONE  FIRST 
C  THEN  THE  MIDDLE  TERMS 

DO  DO  1*1, COUNT 

C  DETERMINING  FIRST  AND  LAST  TERMS  OF  UNLINKED  DATA  FIELD 

IF  ( D I GRAY (1,1)  EO.  ID  .AND.  DIGRAY'  I- 1  , 1  )  EQ .  S35S ) 

+  THEN 

AREAS*' DIGRAYi  I , 3  > *< 0 1  GRAY i  I+NCT-1  , 2 1-0 1  GRAY'! 1*1,0)11 
AREAS*  s:  D I  GRAY  r  I+NCT-1  ,  3  1  S'  D I  GRAY  <  I*NCT-D  ,  0  l-OIGRAY  <  1,2)  O 
GO  TO  DO 
END  IF 

C  DETERMINING  THE  MIDDLE  TERMS 

IF  i:  0 1  GRAY '  1,1)  .  E'J  .  ID  AND  DI  GRAY  (  1*1  ,  1)  .  E'3  .  33?9  > 

THEN 

GO  TO  20 
END  IF 

IF  <: 01  GRAY'!  1,1)  EQ  10)  THEN 

AREAS*  CD  I  GRAY  <!  1 , 3  !>  *  ( 0 1  GRAY  (  1-1  ,2)-DIGRAYf.  1  +  1,2).') 
AREA4=AREA4  +  AREAS 
END  IF 
CONTINUE 

AREAT* '  i!  AR.EA4  ♦  AREAS  +  AREAS  )  / '.  -2 . 0  )  )*1 01  354 1  2 
OPEN'.  UNIT*2 ,  FILE*  ’  PRN '  ,  STATUS*  ’  OLD  •  ) 

WRITEC2, 101 )  (ID) 

WRITE'!  2 ,  * )  (’YOUR  AREA  IN  CUBIC  METERS  IS:  ’) 

WRITE (2, 102) AREAT 
FORMAT  <!  FI  1.2) 

RETURN 
END 

SUBROUTINE  STATS 
IMPLICIT  REAL  (A-2) 

C  THIS  PROGRAM  WILL  COMPUTE  THE  MAXIMUM,  MINIMUM,  AVERAGE,  RANGE,  AND 
C  STANOARO  DEVIATION  FOR  A  GIVEN  DATA  FILE. 

CHARACTER* 12  8EG I N 

PRINT*,  'WELCOME  TO  STATS.  THIS  PROGRAM  WILL  DETERMINE  THE  ' 
PRINT*,  'MAXIMUM,  MINIMUM,  AVERAGE,  RANGE,  ANO  STANDARD  DEVIATION- 
PRINT*,  'FOR  A  GIVEN  FILE  OF  DATA.  YOUR  DATA  FILE  MUST  BE  ONE  ' 
PRINT*,  'STRING  OF  NUMBERS  IF  YOU  ARE  READY,  P LEASE  TYPE  IN  ' 

PRINT*,  'YOUR  FILE  NAME  PLEASE  LIMIT  IT  TO  12  LETTERS. 

REA0(5, 1001BEGIN 

0PEN(UNIT*1 ,  FILE=6EGIN.  STATUS* ' OLD ' ) 

100  FORMAT  <!  A1  2) 

AVG*0  0 

e  i<a*o.o 

SMALL* 1 00000 . 0 
COUNT *0  0 

1  0  REAOl'  1  ,  *  ,  ERR*SOO  '  DATA 


20 


1  02 
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IF  ■  DATA  GT  BIG)  BIG  a  DATA 
IF  i  DATA  LT.  SMALL)  SMALL  =■  DATA 
COUNT  a  COUNT  +  1 
AVG  =  AVG  +  DATA 
GOTO  10 

PRINT*,  'MAXIMUM  VALUE  IS  '.BIG 

PRINT*,  'MINIMUM  VALUE  IS  '.SMALL 

AMEAN  *  AVG  /  COUNT 

PRINT*,  'THE  AVERAGE  IS  '.AMEAN 

RANGE  *  BIG  -  SMALL 

PRINT*.  'THE  RANGE  IS  '.RANGE 

REUINOC 1 ) 

OEVIA  a  0.0 

READ*.  1  ,*,  ERR  »  SOI)  DATA 
SUM  »  f  i  DATA  -  AMEAN )*.*2> 

OEVIA  *  OEVIA  +  SUM 
GOTO  20 

FTOVIA  =  SORT <: OEVIA  /  COUNT) 

PRINT*,  'STANDARD  DEVIATION  IS  ‘.STDVIA 
RETURN 


